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SECTION  I 
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INTRODUCTION 

The  SURVIVE  computer  model  may  be  used  to  evaluate  the  probability  of  sur- 
vival of  a single  penetrator  flying  a specified  flight  path  in  an  environment 
defended  by  surface-to-air  missile  systems  (SAMS).  Options  of  scenario,  coor- 
dinate system,  SAM  firing  doctrine,  and  target  location  give  the  program  flexi- 
bility in  the  types  of  problems  it  is  able  to  handle.  At  one  extreme,  SURVIVE 
can  evaluate  the  survival  probability  of  a weapon  launched  against  a fixed 
target  defended  by  a single  SAM;  at  the  other,  it  can  evaluate  the  expected 
survival  probability  of  a penetrator  and  the  weapon  it  launches  through  a cor- 
ridor defended  by  many  SAMs  of  up  to  10  types. 

The  defense  environment  in  SURVIVE  may  be  specified  in  one  of  two  ways. 
First,  the  exact  locations  of  all  the  SAMs  can  be  specified.  The  model  then 
determines  the  nunber  of  missiles  each  SAM  is  able  to  fire  at  the  penetrator 
during  the  time  it  is  in  coverage  and,  subsequently,  the  penetrator  survival 
probability  based  on  the  single  shot  kill  probability  of  each  SAM.  In  addition 
to  this  specific  approach,  the  model  can  calculate  an  expected  survival  probabi- 
lity by  generating  a representative  sample  of  locations  of  the  SAM  sites  with 
respect  to  the  penetrator  flight  path.  The  model  then  determines  the  number  of 
missiles  each  site  can  fire  at  the  penetrator  and  calculates  an  average  for  all 
locations  where  intercepts  occur.  The  expected  number  of  encounters  that  the 
penetrator  will  have  with  each  SAM  and  the  expected  value  of  the  survival  proba- 
bility are  then  calculated. 

The  SURVIVE  model  incorporates  the  following  aspects  of  the  survivability 
problem: 


i 
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SCENARIO 


Geometry 
SAM  placement 

Masking  of  the  penetrator  by  terrain 
PENETRATOR  AND  WEAPON 

Flight  profiles  (3-dimensional  time  dependent) 

Radar  and  IR  signatures 
Electronic  countermeasures 

SAM  CHARACTERISTICS 
Salvos  per  site 
Missiles  per  salvo 

Timing  criteria  for  tracking,  firing,  and  reloading 

Lethal  envelope 

Missile  flyout  profile 

Kill  probabilities,  single  shot 

Radar  antenna  height 

Radar  ground  clutter  angle 

Radar  sensitivity  or  maximum  range 

Maximum  radar  elevation 

Missile  guidance 

Geometric  launch  restrictions 

ECM  effectiveness 

Probability  of  engagement 

The  SURVIVE  model  is  described  in  section  2;  the  inputs  to  the  model  and  its 

output  in  section  3.  Appendix  A is  a FORTRAN  listing  of  the  program  and  appendix 
B is  a sample  problem. 
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SECTION  2 


MODEL  DESCRIPTION 

A.  GENERAL  DESCRIPTION  OF  THE  MODEL.  Two  methods  are  applicable  to  calculating 
the  probability  of  survival  of  a penetrator  flying  over  an  area  defended  by  SAM 
sites:  the  deterministic  method  and  the  probabilistic  method.  In  a purely 
deterministic  method,  all  the  aspects  of  the  problem  that  are  considered  would 
be  simulated  by  the  model  to  determine  If  the  penetrator  survived  under  a speci- 
fic set  of  circumstances  or  not.  Many  cases  would  have  to  be  run  to  obtain  a 
statistically  valid  sample  from  which  the  probability  of  survival  could  be  deter- 
mined. In  a purely  probabilistic  method,  the  probability  functions  for  the 
various  aspects  of  the  problem  would  be  determined  and  combined  to  obtain  the 
probability  of  survival.  This  tends  to  obscure  the  effects  that  specific  aspects 
of  the  problem  have  on  the  results,  but  simplifies  modeling  of  the  problem.  The 
SURVIVE  program  Is  a hybrid  model.  Many  of  the  Interactions  between  the  SAMs 
and  the  penetrator  are  simulated  deterministically  and  others  are  handled  proba- 
bilistically. For  example,  the  number  of  missiles  that  a SAM  can  fire  at  a pen- 
etrator Is  calculated  deterministically,  but  whether  the  site  will  actually  fire 
Is  controlled  by  a probability  Input  to  the  model.  Two  approaches  are  used  for 
locating  the  flight  profiles  relative  to  the  SAM  sites:  the  specific  approach 
and  the  expected  value  approach. 

In  the  specific  approach,  the  exact  locations  of  the  SAM  sites  are  specified 
In  a two-dimensional  coordinate  system.  A flight  profile  for  the  penetrator  Is 
defined  and  the  number  of  missiles  each  site  can  fire  at  the  penetrator  Is  cal- 
culated. This  value  Is  combined  with  the  single-shot  kill  probability  to  obtain 
the  survival  probability. 

Under  battlefield  conditions  It  is  most  difficult  to  pinpoint  the  exact 
locations  of  SAM  sites.  Intelligence  concerning  their  locations  Is  usually  out- 
dated and  scanty.  Such  uncertainty  In  the  locations  of  the  SAM  sites  relative 
to  the  penetrator  flight  path  requires  a number  of  cases  to  be  run  and  averaged 
to  give  a representative  survival  probability.  This  would  prove  to  be  both 
tedious  and  time  consuming  if  done  manually.  The  expected  value  approach  gene- 
rates a representative  sample  of  geometries  between  the  penetrator  flight  path 
and  the  SAM  sites.  The  number  of  missiles  fired  by  the  SAM  sites  Is  determined 
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for  each  geometry.  The  model_then  computes  the  survival  probabilities  for  each 
of  these  geometries  and  averages  them  to  provide  a representative  survival  proba- 
bility. 

B.  DETAILED  CAPABILITIES  OF  THE  MODEL.  This  section  discusses  the  various 
aspects  of  the  SURVIVE  model  in  detail.  The  program  uses  identifying  labels  on 
the  input  cards  to  associate  input  parameters  with  their  function  in  the  model. 
The  identifier  may  be  associated  with  a single  input  parameter,  in  which  case 
the  identifying  label  is  the  parameter  name.  The  identifier  may  also  be  asso- 
ciated with  a functionally  related  group  of  input  parameters.  In  this  case  the 
identifier  is  a mnemonic  device  to  aid  the  user  in  formatting  the  input  para- 
meters. The  input  parameters  and  identifiers  discussed  in  section  3 are  given 
in  parentheses  throughout  this  section  when  the  concepts  to  which  they  relate 
are  discussed.  Input  parameters  are  also  cross-referenced  to  the  identifier 
with  which  they  are  associated. 

As  mentioned  previously,  the  model  may  be  used  in  its  specific  mode  when  SAM 
placement  with  respect  to  the  penetrator  flight  profile  is  known,  or  it  can  be 
used  in  the  expected  value  mode  when  less  is  known  about  the  placement  of  the 
SAMs  (AVERAGE).  The  expected  value  approach  is  more  useful  in  developing  system 
evaluation  criteria,  since  the  resulting  survival  probability  is  an  average 
based  on  a representative  number  of  relative  placements  of  the  SAMs  with  respect 
to  the  penetration  flight  path.  The  specific  mode  is  most  effective  in  situa- 
tions where  operational  information  would  give  some  idea  of  SAM  locations,  thus 
indicating  a preferred  flight  path  to  avoid  those  SAMs. 

1.  Scenario.  Rectangular  or  polar  coordinates  may  be  selected  (AREA). 
Rectangular  coordinates  are  most  useful  in  describing  the  penetration  of  an 
area  or  a corridor  where  there  is  a preferred  orientation  of  the  flight  path. 
Polar  coordinates  are  more  useful  in  describing  an  attack  on  a point  target 
where  the  penetrator  may  come  from  any  direction. 

Figure  1 shows  the  rectangular  coordinate  system  used  by  the  model.  The 
SAM  sites  are  located  in  X and  Y (XSITE,  YSITE)  when  the  model  is  used  in  its 
specific  mode.  When  the  model  is  run  in  the  expected  value  mode  the  density  of 
the  SAM  sites  must  be  specified.  This  is  done  by  defining  a corridor  of  speci- 
fied width  (CORWDTH).  The  area  defended  by  SAMs  is  then  defined  by  locating 


Figure!.  Rectangular  Coordinates 


the  SAMs  within  the  corridor  as  a function  of  depth  (XSITE  - see  SITE).  The  cor- 
ridor is  centered  about  the  X axis.  The  SAM  is  able  to  assume  any  Y position 
across  the  corridor  for  a given  X.  If  more  than  one  SAM  site  (NSITE  - see  SITE) 
of  the  same  type  (NTYPE  - see  SITE)  is  specified  at  that  X position,  the  sites 
are  assumed  to  be  evenly  spaced  across  the  corridor  to  assure  maximum  coverage 
of  the  penetrator  path.  Two  boundary  conditions  may  be  selected  (SYMETRY).  One 
boundary  condition  describes  a corridor  that  is  a uniform  slice  out  of  an  infinite 
area  with  a homogeneous  defense  structure,  i.e.,  the  defenses  in  the  corridor  are 
repeatedly  duplicated  on  each  side  of  the  corridor.  The  other  describes  an  iso- 
lated corridor  with  no  SAMs  outside  the  corridor. 

Figure  2 shows  the  Polar  Coordinate  System.  The  origin  coincides  with 
that  of  the  rectangular  system.  Angles  are  measured  counter-clockwise  from  the 
X axis.  Site  coordinates  are  specified  by  R and  e (YSITE,  XSITE).  When  the 
model  is  used  in  the  expected  value  mode,  a corridor  must  be  defined  to  specify 
the  density  of  the  SAM  sites.  In  keeping  with  the  polar  coordinates,  the  corri- 
dor is  a sector  originating  at  the  origin  of  the  polar  coordinate  system.  The 
angular  width  of  the  corridor  (CORWDTH)  must  be  specified.  SAM  sites  are  then 
positioned  in  the  corridor  by  specifying  R.  The  angular  position  of  the  center 
of  the  corridor  (TARGETY)  must  be  specified  when  an  isolated  corridor  is  con- 
sidered. This  is  necessary  to  define  the  orientation  between  the  corridor  and 
the  penetrator  flight  path.  The  treatment  of  boundary  conditions  is  similar  to 
the  rectangular  coordinate  case. 

Terrain  masking  of  the  SAM  sensor  (TERRAIN)  is  handled  by  the  model 
through  a discrete  series  of  paired  forward  (in  the  direction  of  decreasing  X) 
and  rearward  masking  angles  (ELF,  ELR  - see  TERRAIN),  and  their  associated  proba- 
bility of  occurrence  (PROB  - see  TERRAIN).  A constant  is  added  to  the  terrain 
angles  to  account  for  the  effect  of  ground  clutter  degrading  sensor  performance 
at  sensor  elevations  just  above  the  terrain  (CLUTTER).  The  differentiation 
between  forward  and  rearward  masking  angles  is  necessary  because  the  SAM  sites 
may  be  backed  up  against  a hill  or  forest  to  protect  them  from  aircraft  coming 
from  the  rear  (figure  3).  Each  pair  of  angles  defines  minimum  elevation  angles 
for  the  sensor  below  which  the  sensor  is  ineffective.  A series  of  masking  angles 
and  their  probability  of  occurrence  may  be  specified  because  in  reality  the  ter- 
rain masking  angles  will  not  be  the  same  for  all  azimuths  in  the  forward  and 
rearward  masking  angle  areas.  The  masking  angles  are  specified  as  a function  of 
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the  sensor  height  above  ground  level  (ANTH  - see  TERRAIN),  as  well  as  terrain 
features  (TFRAC  - see  TERRAIN).  Up  to  10  different  series  may  be  used  to 
describe  different  sensor  heights  and  terrain  features.  The  problem  Is  run  for 
each  pair  of  masking  angles  In  a series.  The  results  are  weighed  by  the  proba- 
bility of  occurrence  to  get  an  average  for  that  series. 

Rectangular  coordinates  are  most  useful  In  describing  scenarios  for 
enroute  attrition  over  a large  area.  Polar  coordinates,  when  used  In  the  expect- 
ed value  mode,  are  convenient  for  certain  types  of  terminal  problems.  For 
example,  consider  a target  defended  by  several  SAMs  that  Is  being  attacked  by  a 
penetrator  that  may  approach  from  any  direction.  This  scenario  Is  easily  de- 
scribed In  the  polar  coordinate  system.  The  target  Is  located  at  the  origin  and 
the  SAMs  at  some  radial  distance  from  It.  The  penetrator  profile  Is  oriented  to 
attack  the  target  at  the  origin  of  the  rectangular  coordinate  system.  The  model 
then  simulates  the  different  directions  of  attack  by  changing  the  positions  of 
the  SAMs  around  the  target  rather  than  by  actually  changing  the  penetrator  pro- 
file. The  probability  of  survival  is  thus  averaged  over  all  directions  of 
attack.  Two  orientations  may  be  selected  for  the  forward  direction  of  the  SAM 
sites  (AREA).  Selecting  the  sites  to  always  Face  radially  outward  from  the 
origin  simulates  a penetrator  approaching  from  any  direction.  Selecting  the 
sites  to  always  face  In  a negative  X direction  of  the  rectangular  coordinate 
system  simulates  an  unknown  angular  location  of  the  sites  around  the  target. 

Once  the  set  of  SAM  sites  (SITE)  has  been  specified,  two  ways  are  avail- 

able for  changing  the  relative  number  of  SAMs  In  the  set.  The  relative  number 
of  each  SAM  type  may  be  altered  by  a multiplicative  factor  (OF).  Additionally, 

It  Is  often  of  Interest  to  look  at  different  multiples  (overall  defense  levels) 
of  the  resulting  set  of  SAM  sites.  The  total  number  of  all  SAMs  specified  for 
the  defense  may  be  altered  by  a multiplicative  factor  (DLEV).  A series  of 
values  may  be  specified  and  the  probability  of  survival  will  be  calculated  for 
each  value. 

2.  Penetrator  and  Weapon.  Two  separate  flight  paths  may  be  handled  by  SUR- 
VIVE at  the  same  time;  those  of  a penetrator  (XYZT)  and  the  weapon  It  launches 

(XYZTW).  The  flight  paths  specify  the  time  dependent  position  of  the  penetrator 
and  weapon  In  three  dimensions.  X and  Y are  In  the  rectangular  coordinate  system 
described  previously,  and  the  third  coordinate  Is  altitude  above  the  X-Y  plane. 


It  is  convenient  to  set  up  the  problem  such  that  the  forward  edge  of  the  battle 
area  (FEBA)  is  located  along  the  Y axis  with  the  area  defended  by  the  SAMs  to 
the  right  of  the  Y axis.  The  SAM  sites  are  assumed  by  the  model  to  be  set  up 
facing  in  the  negative  X direction.  Flight  profiles  would  normally  originate 
to  the  left  of  the  FEBA  and  ingress  to  the  defended  area.  (The  profile  could 
alternately  be  oriented  to  attack  from  behind  the  SAM  sites.) 

It  is  convenient  when  the  model  is  used  in  the  expected  value  mode  to 
specify  the  flight  paths  relative  to  an  imaginary  target  located  at  X and  Y co- 
ordinates of  zero.  The  flight  paths  are  then  shifted  by  an  increment  in  X 
(XSTART)  to  position  the  target  (and,  hence,  the  flight  path)  at  the  desired 
location  in  the  corridor.  The  resulting  probabilities  of  survival  for  any  given 
target  location  are  highly  dependent  on  the  relative  positions  of  the  flight 
paths  to  the  SAM  locations.  Sometimes  it  is  desirable  to  have  an  average  proba- 
bility of  survival  for  a target  situated  over  a range  of  locations.  The  model 
will  generate  a series  of  flight  paths  starting  with  XSTART  and  incrementing  it 
(DXSTART)  to  obtain  a specific  number  of  paths  (NXSTART).  The  program  then 
averages  the  probabilities  of  survival  for  all  target  locations. 

SURVIVE  can  simulate  SAM  radar  and  IR  sensor  performance.  The  model  has 
provisions  for  specifying  radar  and  IR  signatures  for  both  the  penetrator  and 
weapon  CSIGNATURE).  The  signatures  are  specified  as  a function  of  aspect  angle 
around  the  penetrator  or  weapon  (figure  4).  The  radar  signatures  are  specified 
as  the  apparent  radar  cross-section  area.  The  model  determines  the  signal 
strength  at  the  sensor  as  the  apparent  radar  cross-section  area  divided  by  the 
distance  between  the  target  and  the  SAM  site  raised  to  the  fourth  power.  The  IR 
signatures  are  specified  as  radiated  power.  The  signal  strength  at  the  sensor 
is  radiated  power  divided  by  the  distance  between  the  target  and  the  SAM  site 
squared.  Minimum  signal  strengths  for  sensor  tracking  (RTRK  - see  SAM)  and  for 
sensor  lock-on  (RLOCK  - see  SAM)  must  be  specified. 

Electronic  countermeasures  are  not  directly  modeled  in  SURVIVE,  but 
these  effects  are  simulated  by  a degradation  factor  (ECME  - see  SITE)  that  re- 
duces the  effectiveness  of  the  SAM  sites  by  reducing  the  single  shot  kill  proba- 
bility. ECM  degradation  of  the  SAM  sites  may  be  applied  selectively  to  different 
portions  of  the  defended  area  by  specifying  intervals  in  X or  R (XECM).  All  SAM 
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sites  with  X or  R coordinates  within  those  intervals  will  be  treated  by  the  model 
as  being  subjected  to  ECM  degradation. 


3.  SAM  Characteristics.  Ten  different  types  of  SAM  sites  may  be  modeled  in 
SURVIVE  (SAM).  The  characteristics  of  each  SAM  site  type  (number  of  missiles  per 
salvo,  number  of  salvos,  etc.)  are  specified  by  input.  Each  type  may  be  differ- 
entiated according  to  its  firing  capabilities,  missile  performance,  sensor  track- 
ing, and  missile  launch  characteristics.  The  firing  capability  of  a SAM  is 
described  by  the  number  of  missiles  it  fires  per  salvo  (NSS  - see  SAM),  the  time 
between  firing  each  missile  in  a salvo  (TISH  - see  SAM),  the  nunber  of  salvos 
that  can  be  fired  before  reloading  (NS  - see  SAM),  the  minimum  time  between  suc- 
cessive salvos  (TINTER  - see  SAM),  the  missile  launcher  reloading  time  (TRELOAD  - 
see  SAM),  and  the  maximum  azimuth  angle  at  which  missiles  may  be  launched,  mea- 
sured from  the  forward  direction  of  the  site  (AZMAX  - see  SAM).  Missiles  are 
assumed  to  fly  straight  line  intercepts  and  flight  times  are  specified  by  a time 
vs  distance  function  for  each  SAM  type  (MISLXT).  Missile  intercept  envelopes 
are  described  by  a minimum  intercept  altitude  (ALTMIN  - see  SAM),  as  well  as 
dead  zone  (FUSE)  and  maximum  lethal  range  (RN6)  as  a function  of  elevation  angle 
(figure  5).  Single  shot  kill  probabilities  (PKSS)  are  used  to  describe  the 
average  lethality  oP  a single  missile  against  the  penetrator  and  the  weapon,  but 
need  not  be  the  same  for  both. 

The  geometric  characteristics  of  the  tracking  system  are  described  by  the 
height  above  ground  level  of  the  antenna  or  other  sensor  (HRAD  - see  SAM)  and  its 
maximum  elevation  angle  (ELMAX  - see  SAM).  Several  restrictions  may  be  placed  on 
the  sensor's  performance.  Maximum  range  constraints  may  be  used  for  specifying 
initiation  of  tracking  (RADTRK  - see  SAM)  and  guidance  system  lock-on  (RLOCK  - 
see  SAM),  or  minimum  signal  strengths  for  tracking  (RTRK  - see  SAM)  and  lock-on 
(RLOCK  - see  SAM)  can  be  used  in  conj'unction  with  penetrator  and  weapon  signa- 
tures to  more  accurately  predict  sensor  performance.  Restrictions  may  also  be 
placed  on  the  aspect  angle  between  the  line  of  sight  vector  from  the  SAM  site  to 
the  penetrator  or  weapon  and  the  velocity  vector  of  the  penetrator  or  weapon  at 
the  time  of  missile  launch  (figure  6).  The  velocity  vector  is  tangent  to  the 
flight  path.  The  aspect  angle  is  zero  when  the  penetrator  or  weapon  is  flying 
directly  at  the  SAM  site,  and  180  degrees  when  flying  directly  away  from  it.  The 
minimum  and  maximum  aspect  angles  for  firing  (ASPMIN,  ASPMAX  - see  SAM)  allow  the 
model  to  simulate  systems  that  cannot  lock  on  to  targets  with  negative  Doppler 
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shifts,  i.e.,  moving  away  from  the  SAM  site.  This  feature  can  also  be  used  to 
simulate  an  IR  system  in  the  absence  of  signature  data  because  in  most  cases  an 
IR  sensor  cannot  lock  on  to  a penetrator  when  it  is  headed  toward  the  sensor, 
i.e.,  for  small  aspect  angles,  since  it  cannot  see  the  heat  source. 

A look-shoot-look  firing  doctrine  is  employed  for  the  SAMs.  The  SAM  must 

be  able  to  track  the  penetrator  or  weapon  for  a specified  time  (TINIT  - see  SAM) 

before  it  can  fire  a salvo.  It  must  wait  the  same  amount  of  time  after  the  last 

missile  of  the  salvo  "intercepts"  to  evaluate  the  effects  of  the  salvo  and  update 

its  tracking  information  before  it  can  fire  again.  This  interval  may  not  be  less 
than  the  minimum  time  between  salvos  or  the  reload  time,  depending  on  the  number 
of  salvos  the  site  has  already  fired.  The  site  will  fire  as  long  as  it  is  able 
to  track  and  intercept  the  penetrator. 

The  model  allows  two  types  of  missile  guidance  systems  (IR  - see  SAM). 

The  first  is  a guidance  system  that  is  self-homing,  such  as  an  IR  seeker;  the 
second  must  have  the  target  in  sensor  coverage  from  the  site  to  guide  the  missile 
to  intercept. 

It  is  possible  that  SAMs  would  not  fire  as  readily  at  egressing  penetra- 
tors,  preferring  to  save  their  missiles  for  ingressing  penetrators  that  might  be 
more  of  a threat.  The  model  can  allow  for  this  (EGRESS)  by  specifying  a time 
(TEGRESS  - see  EGRESS)  corresponding  to  the  point  on  the  penetrator  flight  path 
after  which  the  SAMs  firing  will  be  less  frequent,  and  a multiplicative  factor 
(FEGRESS  - see  EGRESS)  that  will  degrade  the  performance  of  the  SAM  sites  after 
that  time. 

Two  operational  characteristics  of  the  SAMs  are  site  specific.  When  the 
locations  and  numbers  of  the  sites  are  specified  (SITE),  the  probability  that 
the  sites  at  each  location  will  fire  at  a penetrator  (PUP  - see  SITE)  and  the 
terrain  identifying  factor  (TERFR  - see  SITE)  for  those  sites  must  be  specified. 

C.  CALCULATION  OF  AVERAGE  NUMBER  OF  MISSILES  FIRED  AND  SITE  LETHAL  WIDTH.  The 
model  simulates  encounters  between  the  SAMs  and  the  penetrator  and  its  weapon  as 
a function  of  time  to  determine  when  missiles  may  be  fired.  Start  and  stop  times 
corresponding  to  the  part  of  the  penetrator  profile  (TSTART,  TSTOP)  and  weapon 
profile  (TSTARTW,  TSTOPW)  to  be  considered  in  the  survivability  calculations 
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must  be  specified.  The  model  will  derive  the  simulation  time  interval  limits 
(TSTARTP,  TSTOPP)  from  the  start  and  stop  times  of  the  penetrator  and  weapon.  If 
a value  of  TSTARTP  or  TSTOPP  is  specified  by  input  to  the  model,  it  will  super- 
sede the  derived  value.  The  model  calculates  both  the  average  number  of  missiles 
fired  and  site  lethal  width  when  used  in  the  expected  value  mode.  When  used  in 
the  specific  mode,  it  calculates  the  number  of  missiles  fired  only  from  the  spe- 
cific site  locations  and  does  not  vary  the  geometry  between  the  sites  and  the 
penetrator  path.  The  site  lethal  width  is  not  calculated. 

In  calculating  the  lethal  width  of  a SAM  site  and  average  number  of  mis- 
siles fired,  the  model  considers  only  one  SAM  site  at  a time.  Using  the  maximum 
range  of  the  missile,  the  model  determines  the  time  intervals  along  the  flight 
paths  that  the  penetrator  and  weapon  are  in  coverage  of  the  site.  It  then  finds 
the  range  in  the  Y direction  from  the  site  for  which  intercepts  can  occur.  This 
range  is  then  subdivided  to  give  a representative  number  of  different  geometries 
by  offsetting  the  flight  path  from  the  site  in  the  Y direction.  The  model  maps 
the  flight  paths  and  SAM  sites  into  a spherical  earth  coordinate  system  to  pro- 
vide a more  accurate  representation  of  the  problem  geometry  by  accounting  for  the 
curvature  of  the  earth.  This  mapping  is  accomplished  by  translating  linear  di- 
mensions in  the  X-Y  coordinates  into  arc  lengths  at  the  corresponding  altitudes 
above  the  surface  of  the  spherical  earth.  This  provides  a representation  of 
horizon  effects  on  the  geometry  between  the  SAM  site  and  the  flight  paths.  The 
model  then  examines  the  in-coverage  time  intervals  using  a small  time  step  to 
determine  the  number  of  missiles  the  site  can  fire  at  the  profiles  for  each  off- 
set. The  lethal  width  of  the  site  is  calculated  by  multiplying  the  number  of 
offsets  with  at  least  one  missile  firing  by  the  offset  subdivision  interval.  The 
number  of  missiles  fired  by  the  site  are  averaged  over  all  offsets  for  which 
intercepts  are  possible.  The  model  keeps  track  separately  of  the  average  number 
of  missiles  fired  and  site  lethal  width  for  the  penetrator  prior  to  weapons  re- 
lease, for  the  complete  flight  path,  and  for  the  weapon  during  its  flight.  The 
sites  will  continue  to  fire  at  the  penetrator  after  the  penetrator  launches  its 
weapon  until  a missile  may  be  fired  at  the  weapon.  At  this  point,  the  site  will 
guide  any  previously  fired  missiles  to  the  penetrator  Cif  the  missiles  are  not 
self-homing)  and  starts  tracking  and  firing  at  the  weapon.  The  site  will  fire  at 
the  weapon  throughout  its  flight.  When  the  SAM  site  is  no  longer  able  to  fire  at 
the  weapon,  it  will  re-acquire  the  penetrator  if  it  is  in  coverage  and  continue 
firing  at  it. 
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The  calculation  of  the  average  number  of  missiles  fired  is  quite  time 
consuming,  and  the  model  makes  several  provisions  for  decreasing  the  calculation 
time  at  the  expense  of  precision.  One  is  a parameter  for  increasing  the  time 
step  for  examining  the  in-coverage  intervals  (FASTRUN).  The  value  multiplies  the 
normal  timestep.  A value  of  three  will  cut  running  time  by  about  50%  and  only 
changes  the  calculated  value  of  survival  probability  by  about  3%. 

The  other  provision  eliminates  the  need  for  calculating  the  average  niro- 
ber  of  shots  and  lethal  width  for  similar  sites  that  cover  identical  portions  of 
the  flight  profiles  when  the  model  is  used  in  its  expected  value  mode.  On  many 
types  of  profiles,  a significant  portion  of  the  flight  path  will  appear  identical 
to  any  SAM  site  of  a given  type  that  can  shoot  only  at  that  portion  of  the  path. 
The  model  has  an  input  parameter  for  each  flight  path  that  specifies  the  largest 
value  of  the  X coordinate  for  this  portion  of  its  flight  profile  (XYZT,  XYZTW). 

For  example,  an  aircraft  flight  profile  that  ingresses  at  constant  speed  and  alti- 
tude performs  a maneuver  and  then  egresses  parallel  to  its  ingress  path  at  a 
constant  speed  and  altitude  would  appear  identical  to  a SAM  site  at  any  location 
up  to  the  area  in  which  the  maneuver  is  performed.  In  figure  7 the  flight  path 
would  appear  identical  to  a SAM  site  located  at  a constant  Y and  at  any  X up  to 
the  point  where  it  could  intev'act  with  the  penetrator  during  its  maneuver.  It  is 
obvious  that  the  average  number  of  missiles  fired  and  site  lethal  width  for  SAMs 
A,  B,  and  C will  be  the  same,  so  they  need  only  be  calculated  once. 

0.  CALCULATION  OF  SURVIVAL  PROBABILITY.  When  the  model  has  calculated  the  aver- 
age number  of  missiles  fired  by  all  SAM  sites  and  their  lethal  widths,  all  the 
information  necessary  for  calculating  the  survival  probability  is  available.  The 
probability  of  survival  is  calculated  for  each  value  of  the  series  of  defense 
levels  (DLEV). 

The  probability  of  a single  penetrator  surviving  the  attack  of  type  j 
SAMs  at  a site  location  i for  defense  level  m is: 
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Figure  7.  Portion  of  Plight  Path  That  Appears  Similar  to  SAMs  at  Different  Locations 
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NSITE^. 

number  of  SAM  sites  at  location  i 

Si 

number  of  missiles  fired  by  a j type  SAM  at  location  i 

fraction  of  total  number  of  j type  SAMs  operative 

PUP^ 

probability  that  a j type  SAM  site  at  i will  fire  at 
the  penetrator,  given  the  opportunity 

defense  level  fraction 

PEi 

expected  number  of  encounters  of  penetrator  with  site 

where 
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specific  placement  of  SAMs 
relative  to  target  path 
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expected  value  encounter  with 
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expected  value  encounter  with 

isolated  corridor  boundary 

condition 

lethal  width  of  SAM  at  i 
corridor  width  (CORWDTH) 


*The  factor  (L.*2W)'^  is  the  expected  number  of  encounters  with  sites  outside  the 
boundaries  of  the  corridor  and  must  be  subtracted  when  the  isolated  corridor 
boundary  is  used. 
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The  probability  of  the  penetrator  surviving  all  the  SAM  locations  for  defense 
level  m is; 
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SECTION  3 


INPUT  AND  OUTPUT 

A.  INPUT.  Input  to  the  SURVIVE  program  is  based  on  cards  containing  an  identi- 
fier and,  generally,  data  associated  with  the  identifier.  The  identifier  may  be 
the  name  of  a variable  in  the  model.  The  first  ten  spaces  on  the  card  are  used 
for  the  identifier,  left  justified,  and  the  second  ten  spaces  are  a floating  point 
field  for  data.  The  remainder  of  the  card  may  be  used  for  comments.  In  many 
instances,  this  format  is  not  adequate  to  allow  the  specification  of  all  data 
associated  with  the  identifier.  In  these  cases,  the  identifier  card  is  itrme- 
diately  followed  by  additional  specially  formatted  data  cards.  The  order  of  the 
identifier  cards  within  a data  deck  is  not  important,  although  cards  associated 
with  an  identifier  card  must  be  in  proper  order.  A specific  case  is  set  up  by 
defining  all  the  necessary  information  for  the  model  with  the  identifier  cards 
and  their  associated  cards.  The  case  is  terminated  by  the  identifier  "ENDCASE." 
Multiple  cases  may  be  run  by  the  model.  Once  data  are  defined  by  an  identifier 
case,  they  may  be  changed  by  specifying  new  data  with  another  identifier  card 
with  the  same  identifier  name.  Data  which  does  not  change  from  case  to  case  need 
only  be  input  once.  Additional  cases  may  be  defined  by  changing  parameters  of 
the  previous  data  set,  thus  simplifying  the  generation  of  multiple  cases.  Program 
execution  is  iitinediately  terminated  by  the  identifier  "ENDJOB."  The  program 
prints  out  all  the  input  data  to  provide  a permanent  record  of  the  parameters 
associated  with  each  run.  Data  appearing  on  cards  associated  with  the  identifier 
card  are  labeled  in  the  output  with  variable  names  or  descriptions.  These  para- 
meter names  and  descriptions  were  given  in  the  previous  section  to  show  their  use 
in  the  model.  In  this  section  they  are  used  in  addition  to  the  identifier  names 
to  help  describe  the  input  to  the  model.  The  model  checks  for  invalid  identifier 
names  and  will  issue  a diagnostic  message  and  terminate  execution  when  one  is 
found.  Input  cards  for  a sample  problem  are  given  in  appendix  B. 

The  following  list  defines  the  identifiers  used  for  input  to  the  model 
along  with  any  other  cards  associated  with  the  identifier.  When  an  identifier 
card  is  used  to  specify  a SAM  type,  the  identifying  number  of  the  SAM  from  1.0 
to  10.0  is  entered  in  the  data  field.  In  some  cases,  a variable  number  of  cards 
may  follow  the  identifier  card,  as  when  defining  a flight  profile.  These  groups 
of  cards  are  terminated  by  an  end  of  record  (EOR)  card.  The  model  reads  data 
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until  the  EOR  card  is  reached,  thus  relieving  the  user  of  the  task  of  counting 
the  cards.  The  units  used  by  the  model  are  kilometers,  seconds,  degrees,  and 
appropriately  derived  units.  Parameter  default  values  are  given  below  when  appli- 
cable. 


Input  Glossary 


AREA 

AVERAGE 

CLUHER 

CORWDTK 

DEBUG 

DEBUGl 

DF 

DLEV 

DXSTART 

EGRESS 


Signifies  rectangular  coordinates  (=  1.),  polar  coordinates 
with  site  facing  to  decreasing  values  of  X (f  Q.)»  or  polar 
coordinates  with  site  facing  radially  outward  from  the  origin 
of  the  coordinate  system  (=  -1.).  Default  value  = 1. 

Specific  mode  (=  0.),  or  expected  value  mode  (=  1.).  Default 
value  = 1. 

Ground  clutter  angle  above  terrain  for  sensors. 

Corridor  width. 

Print  debugging  information  concerning  relative  positions  of 
sites  and  targets  at  each  time  step  (=  1.),  no  information 
printed  (=  0.).  Default  value  = 0. 

Print  debugging  information  concerning  relative  positions  of 
sites  and  targets  for  each  missile  launch  (=  1,),  no  infor- 
mation printed  (=  0.).  Default  value  = 0. 

The  i-th  field  of  the  next  card  contains  the  fraction  of 
type  i SAMs  that  are  to  constitute  the  actual  defense  level. 
Format  (10E8.1). 

Each  following  card  defines  an  overall  defense  level  multi- 
plier for  editing  the  final  output.  Format  (E10.3).  Termi- 
nated by  an  EOR. 

Increment  in  X for  generating  a series  of  delivery  system 
profiles.  Default  value  = 2.5. 

The  following  card  provides  information  on  egress  time  of  the 
penetrator  and  SAM  degradation  factor;  TEGRESS  the  time  after 
which  the  penetrator  is  assumed  to  be  egressing  and  FEGRESS 
the  degradation  factor  applied  to  SAM  performance  after  that 
time.  Format  (2E10.3). 
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ENDCASE 

ENDJOB 

FASTRUN 

FUSE 

MISLXT 

NOSIG 

NXSTART 

PKSS 

RNG 

SAM 


Signifies  the  end  of  input  for  a case. 

Immediately  terminates  execution. 

Multiplies  the  normal  program  time  step  to  decrease  running 
time.  Default  value  = 1. 

Specifies  a SAM  type.  The  following  cards  define  the  dead 
zone  about  the  site  by  paired  elevation-range  points  ordered 
by  increasing  elevation.  Format  (2E10.3).  Terminated  by  an 
EOR. 

Specifies  a SAM  type.  The  following  cards  define  the  missile 
intercept  performance  by  paired  distance-time  points  along 
its  flight  path  ordered  by  increasing  distance.  Format 
(2E10.3).  Terminated  by  an  EOR. 

The  i-th  field  of  the  following  card  indicates  that  for  SAM 
type  i signature  information  and  sensor  sensitivities  (=  0), 
or  maximum  sensor  ranges  (=  1)  will  be  used  for  tracking  and 
lock-on.  Format  (1018).  Default  values  = 1. 

Total  number  of  flight  profiles  to  be  generated  for  this 
case.  Default  value  = 1. 

Selects  the  penetrator  (=  0.)  or  the  weapon  (=  1.),  and  spe- 
cifies that  the  single  shot  kill  probabilities  for  the  type 
i SAM  against  that  target  are  given  in  the  i-th  field  of  the 
following  card.  Format  {10E8.1). 

Specifies  a SAM  type.  The  following  cards  define  the  missile 
maximum  lethal  range  envelope  by  paired  elevation-range 
points  ordered  by  increasing  elevation.  Format  (2E10.3). 
Terminated  by  an  EOR. 

Specifies  a SAM  type.  The  following  three  cards  define 
various  parameters  for  that  type  SAM  site  as  follows: 


Card  1:  Format  (3no,5E10.3) 

field  (1)  NS  number  of  salvos  that  the  site  Is 

able  to  fire  before  reloading. 

(2)  NSS  number  of  missiles  per  salvo. 

(3)  IR  missile  guidance  Independent  of  site 

(=1),  target  must  remain  within 
radar  coverage  of  site  during  mis- 
sile flight  {»  0). 

(4)  HRAD  height  of  sensor  above  ground  level. 

(5)  RTRK  minimum  signal  strength  for  sensor 

tracking. 

(6)  RLOCK  minimum  signal  strength  for  sensor 

lock-on. 

(7}  ELNAX  sensor  maximum  elevation. 


i 


(8)  ALTMIN  minimum  missile  Intercept  altitude 


Card  2:  Format  (8E10.3) 

field  (1)  UNIT  sensor  tracking  time  before  each 

salvo  Is  fired. 


(2)  TISH  time  between  shots  within  a salvo. 


(3)  TINIER  time  between  salvos.  j 

(4)  TRELOAD  site  reload  time.  I 

i 

(5)  ASPMIN  minimum  aspect  angle  of  target  for 

sensor  to  acquire  penetrator  or 
weapon . 

(6)  ASPMAX  maximum  aspect  angle  of  target  for 

sensor  to  acquire  penetrator  or 
weapon. 

(7)  AZMAX  maximum  azimuth  angle  for  firing  a 

missile. 

(8)  ECME  effectiveness  of  site  when  sub- 

jected to  ECM. 
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SIGNATURE 


SITE 


Card  3:  Format  C2E10.31 

field  01  RADTRK  maximum  sensor  tracking  range. 
(2)  RAOLOCK  maximum  sensor  lock-on  range. 


Signifies  signature  data  for  (=  1.)  the  penetrator  for  use  by 
those  sensors  with  IR  = 0,  for  (=  2.)  the  penetrator  for  use 
by  those  sensors  with  IR  = 1,  for  (=  3.)  the  weapon  for  use 
by  those  sensors  with  IR  = 0,  and  for  (=»  4.)  the  weapon  for 
use  by  those  sensors  with  IR  = 1.  The  signature  data  is 
specified  on  the  succeeding  cards  as  a function  of  aspect 
angle  in  order  of  ascending  angles.  The  angle  is  the  first 
value  on  each  card  and  the  signature  second.  Format  (2E10.3). 
Terminated  by  an  EOR. 


The  cards  that  follow  define  the  entire  basic  set  of  SAMs 
available  for  the  case,  their  number,  locations,  and  some 
site  specific  parameters.  Only  one  SAM  type  at  one  location 
may  be  specified  per  card,  but  more  than  one  site  may  be 
specified  at  that  location.  Each  card  defines  the  following 
variables: 


(1)  NTYPE 

(2)  NSITE 

(3)  XSITE 

(4)  YSITE 


(5)  PUP 


type  number  of  SAM. 

number  of  sites. 

X or  0 coordinate  of  sites. 

Y or  R coordinate  of  sites.  It  is 
not  necessary  to  specify  Y or  9 
when  running  in  the  expected  value 
mode  (AVERAGE  =1.). 

probability  that  the  sites  will  fire 
given  an  opportunity. 


(6)  TERFR  terrain  identifying  factor  for  se- 

lecting the  proper  terrain  masking 
angle  distribution  (see  TERRAIN). 
Format  (2110,  4F10.3)  terminated  by  EOR 
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SYMETRY 

TARGETY 

TERRAIN 


TITLE 

TSTART 

TSTARTP 

TSTARTW 

TSTOP 

TSTOPP 


Signifies  Komogeneous  boundary  condition  (=  1.0)  or  isolated 
area  boundary  condition  (=  0.).  Default  value  = 1. 

The  Y value  of  the  ground  target  coordinate  relative  to  the 
delivery  system  profile  when  AREA  = 1,,  or  the  central  angle 
of  an  angular  corridor  when  AREA  = 0.  Only  used  when 
SYMETRY  = 0.  Default  value  = 0. 

Signifies  that  groups  of  discrete  terrain  masking  angle  in- 
formation will  follow.  The  first  card  of  each  of  the  groups 
following  the  TERRAIN  card  specifies  TFRAC,  the  terrain  mask- 
ing angle  identifying  factor,  and  ANTH,  the  sensor  height 
above  ground  level  for  which  the  masking  angles  were  gene- 
rated. Format  (2E10.3).  The  remaining  cards  in  each  group 
specify  the  forward  and  rear  masking  angles,  ELF  and  ELR,  as 
well  as  the  probability  that  they  will  occur,  PROS.  Format 
(3E10.3).  A maximum  of  ten  groups  may  be  specified  with  up 
to  10  paired  angles  in  each  group.  Each  group  is  terminated 
by  an  EOR,  and  an  additional  EOR  must  appear  after  the  last 
group.  The  model  chooses  the  appropriate  group  of  terrain 
masking  angles  to  match  the  terrain  identifying  factor  (TERFR) 
given  on  the  site  location  cards  (SITE),  as  well  as  the 
height  of  the  sensor  for  the  particular  type  of  site  (HRAD) 
as  specified  by  the  site  characteristics  (SAM).  It  is  neces- 
sary to  have  groups  of  terrain  masking  angles  for  all  result- 
ing combinations  of  HRAD  and  TERFR. 

Any  remark  punched  in  the  comment  field  of  this  card  will  be 
used  to  title  the  output. 

Starting  time  of  the  penetrator  flight  profile. 

Simulation  starting  time.  If  not  specified,  the  model  will 
choose  the  smaller  of  TSTART  and  TSTARTW. 

Starting  time  of  the  weapon  flight  profile. 

Stopping  time  of  the  penetrator  flight  profile. 

Simulation  stopping  time.  If  not  specified,  the  model  will 
choose  the  larger  of  TSTOP  and  TSTOPW. 
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TSTOPW  Stopping  time  of  the  weapon  flight  profile. 

XECM  Signifies  that  the  following  cards  will  specify  the  starting 

and  ending  values  of  intervals  in  X or  R over  which  the  SAM 
sites  will  be  degraded  by  ECM,  with  one  interval  per  card. 
Format  (2E10.3).  Terminated  by  an  EOR. 

XSTART  Constant  added  to  the  X coordinate  of  the  penetrator  and 
weapon  flight  profiles  for  locating  it  relative  to  the  SAM 
sites.  Will  be  the  first  value  when  generating  a series  of 
flight  profiles  for  a given  case. 

XYZT  Specifies  the  largest  X coordinate  for  the  portion  of  the 

penetrator  flight  profile  that  will  appear  identical  to 
similar  SAM  sites.  A value  of  zero  indicates  that  all  SAM 
sites  will  be  calculated  individually.  The  following  cards 
specify  penetrator  X,  Y,  Z,  and  corresponding  time  for  each 
point  of  the  flight  profile  in  order  of  increasing  time. 
Format  (4E10.3).  Terminated  by  an  EOR. 

XYZTW  The  same  as  XYZT  except  for  the  weapon  flight  profile. 

B.  OUTPUT.  After  listing  all  data  input  to  the  model,  survival  probabilities 
and  related  information  for  each  case  are  reported.  Figure  8 shows  part  of  the 
survival  probability  results  for  the  sample  problem  in  appendix  B.  Various 
features  of  the  results  will  be  identified  in  figure  8 as  they  are  discussed 
in  the  text. 

The  model  sumnarizes  the  results  for  each  XSTART  generated  for  a given 
case  (A).  For  each  of  the  overall  defense  levels  input  to  the  model,  survival 
probabilities  are  reported  without  any  ECM  degradation  of  the  SAM  sites  (C),  as 
well  as  with  any  ECM  degradation  specified  (D).  Under  each  of  these  headings 
survival  probabilities  as  a function  of  overall  defense  level  are  given  for:  (1) 
the  penetrator  from  the  simulation  start  time  to  the  weapon  start  time  (release) 
(E);  (2)  the  penetrator  from  the  simulation  start  time  to  the  simulation  stop 
time  (F);  and  (3)  the  weapon  from  its  start  time  to  its  stop  time  (G).  In  the 
last  column  (H)  is  given  the  combined  weapon  survival  probability  composed  of 
the  penetrator  survival  probability  with  ECM  up  to  weapon  release  and  the  weapon 
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survival  probability  without  ECM.  The  site  lethal  width  and  average  number  of 
missiles  fired  for  a single  SAM  site  that  are  used  to  calculate  the  survival 
probabilities  are  given  by  SAM  type  and  location  (I). 


After  these  results  have  been  given  for  all  values  of  XSTART  generated 
for  the  case,  the  average  survival  probabilities  for  the  case  are  given  (B). 


APPENDIX  A 


FORTRAN  LISTING  OF  SURVIVE 


I 
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PROGRAM  SURVIVE < INPUT *OUTPUT « T APEl > INPUT ) 


THIS  ROUTINE  SEQUENCES  THE  PROGRAM  FLOW  AND  PRINTS  CASE  AVERAGES 


COMMON  /TRAJ/  NXYZT (2) ,XYZ < lOOOtSt 2) *T ( 1000 *2) •XSTART*TSTART (2) tTS 
1T0P(2) tTINf 10 t2)*NTIN«XECM (20*2) «NXECM«THIN( 10*2) »OXSTART*NXSTARTf 
2XSAMEP*XSAMEW»NSYS0PtTST0PPf TSTARTP,TARGETY*0XYZ0T (lOOOtSf 2) • IS (2> 
COMMON  /PARM/  OF ( 1 0 ) t DEBUG fFLTWTH(2*2)»PLETH(lO)* OEBUGl « T I TLE ( 6 ) • P 
1KSS(10«2>»CORWDTH»CLUTTER«TERRANE(10*3«10)«NTERA(10) *NTER*TFRAC(10 
2 ) * ANTH ( I O ) « ARE A • 0EBU62  « OFSETO •RELEASE • VELPEN , OLE V ( 1 0 ) • NOLE V • A VEP AG 
3E*NITL1<6)«NITL2(6)»SYMETRY»TEGRESS«FEGRESS»NITL3(6)*NINTRPR(10)«N 
AASPvFASTRUN 

COMMON  /AVG/  SUMECM(10,3) •SUMNECM(10»3) 

1 CALL  INPUTS 
DO  2 I«l»30 

2 SUMECMnVsSUMNECMIDsO. 

IF  fNXSTART.LE.O)  NXSTaRT>1 

SVXS«XSTART 

00  3 I>1»NXSTART 

CALL  AVSHOTS 

CALL  PROBS 

XSTART»XSTART»DXSTART 

3 CONTINUE 
XSTARTsSVXS 

PRINT  4,  NXSTART«XSTART 

PRINT  5*  TITLE»NiTLl»NITL2»NITL3 

PRINT  6«  AREA.CORWOTH,XSTART*TSTART,TSTOP 

PRINT  8 

PRINT  7*  fOLEV<L>*(SUMNECM(L«I)«I«l*3)* (SUMECMIL*!) •l3l*3)»SUMECM( 
1L»1)*SUHNECM(L*3)«Ls1*N0LEV) 

PRINT  9 
60  TO  1 
C 

4 FORMAT  C//130<1H*)//22H  AVERAGES  FOR  PREVIOUS. I5*22H  CASES  FIRST 
IXSTART  >«F10«2> 

5 FORMAT  flX«6AI0) 

6 FORMAT  f6H  AREAa*F5.1*l5H  CORRIDOR  WIOTH*F10.2«8H  XSTART3>F10.2/12 
IH  START  TIMES»2FI0.2.12H  STOP  TIMES.2F10.2) 

7 FORMAT  (8F15.3) 

8 FOPMAT  (21X»32HSURVIVAL  PROBABILITY  WITHOUT  ECM*23X«29HSURVIVAL  PR 
lOBABiLITY  WITH  ECM/2X* 13H0EFENCE  LEVEL*? < IX* 14HA/C  TO  RELEASE*3X*1 
22HA/C  complete *5X*10HWEAPON/REL)2X*13HA/C«W  W/0  ECM) 

9 FORMAT  (130(1H*)) 

END 
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SUBROUTINE  AVSHOTS 


THIS  ROUTINE  DETERMINES  OFFSET  RANGES*  CALCULATES  NUMBER  OF  SHOTS 
FIRED*  AND  THE  AVERAGE  NUMBER  OF  SHOTS  ANO  LETHAL  MIOTH 


COMMON  /TRAJ/  NXYZT (2) ,XYZ C 1000*3*2) *T ( 1000*2) *XSTART*TSTART 12) *TS 
1T0P(2)*TIN<10»2)*NTIN»XECM(20*2) *NXECH*THIN( 10*2) *OXSTART*NXSTART* 
2XSAMEP*XSAMEM*NSYS0P*TST0PP*TSTARTP*TARGETY*DXYZ0T (1000*3*2)* IS (2) 
COMMON  /ISAM/  NTYPE(IOO)  *NSITEn00)  *XSITE(100)  »PUP<100)  *ELMIN(  <00* 
12)*NTOTS*SITWOTH(100*3)*AVSHOT(100*3) *TERFR( 100) *SITRAO( 100) *ARSIT 
2(100) 

COMMON  /ASAM/  HRAO (10) »RTRK (10) *ELMAX (10>*TINIT(10) *TINTER ( 1 0 ) *NS ( 
110)  vNSSdO)  •TRELQAO(IO)  *AVVEL(10)  *ASPMIN(10)  *ASPMAX(10)  *AZMAX(10)  * 
2RN6(20*10) •ELR(20*10)*NRN6(10)*FUS(20*10>*ELP(20*10) *NFUS(10) *IR(1 
30)*TISH(10>*ECME(10)*AlTHIN(10)*ALTMAX(10)*SIGTH(20*4)«SIG(20*4)*N 
4SI6(4)  yRLOCKdO)  *XMISL (20*  10)  *TMISL  (20*  10)  *NXMISL ( 10 ) •RAOTRK  (10)  «R 
SAOLOCKdO) 

COMMON  /P ARM/  OF (1 0 ) * OEBUB* FLT WTH ( 2*2)«RLETH(10)*OEBUG1*TITLE(6)*P 
1KSS(10*2)  «CORHOTH»CLUTTER*TERRANE(10*3*10)  *NTERAdO)  *NTER*TFRACdO 
2 ) * ANTH (1 0 ) *AREA • 0EBU62 • OFSETO «RELE ASE  * VELPEN • OLE V (1 0 ) « NOLEV  * AVER AG 
3E*NITL1(6)*NITL2(6)*SYMETRY*TEGRESS*FEGRESS*NITL3(6)*NINTRPR(10)*N 
4ASP«FASTRUN 

COMMON  /O/  AZT *ELT *RN6T* ASPT * AZTF* ELTF *RN6TF • ASPTF * AZTA*ELTA *RNGTA 
1*ASPTA 

DIMENSION  NOONE(IOO)*  NSHT(3)*  AVSH(3)*  N0FSET(3) 

00  1 1*1*100 
1 NOONE(I)90 
0PR-57.29S78 
NS0T«FASTRUN-1. 

IF  (NSOT.LT.O)  NSOTsO 

NOEG«10 

NOOFS*10 

NOTS*50 

TINT-15, 

SVTl-TSTARTP 
SVT2-TST0PP 
00  45  I*1*NT0TS 
TSTARTP«SVTl 
TST0PP*SVT2 

IF  (NOONE(I),GT.O)  60  TO  45 
XSITS«XSITE(I) 

N*NTVPE ( I ) 

SITWOTH(I*l)sSITWOTH(I,2)*SITWOTH(I*3)«0, 

AVSHOT(I*l)sAVSHOT(I*2)*AVSHOT(I*3)*0, 

IF  (OF(N).LE.O,)  go  to  45 
IF  (NTER.LE.O)  GO  TO  4S 
IS(1)*IS(2)«2 
OT*S, 

NSYSOP-O 

IF  (NXYZT(l) •GT,0.A.PKSS(N*1) .GT.O,.A. (IR(N).LE.O.A.NSIG(I) .GT.0,0 
1,IR(N),GT,0,A,NSIG(3),GT,0,0,NINTRPR(N) .GT,0) ) NSYSOP-l 
IF  (NXYZT(2),GT.0.A.PKSS(N«2) .GT^O^.A, ( IR(N) .LE.O.A.NSIGIZ) .GT.0.0 
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1.IR<N).6T.O.A.NSI6(3>.6T.O.O.NINTRPR(N).GT.O) ) NSYS0PsNSYS0P«2 
IF  (NSYSOP,LE*0)  GO  TO  45 
IF  (TSTARTP.GT.-IOOOOO.)  GO  TO  5 
GO  TO  (2»3t4)»  NSYSOP 

2 TSTARTP*TSTART 
60  TO  5 

3 TSTARTP»TSTART(2) 

GO  TO  5 

4 TSTARTP«AMIN1 (TSTART ( I ) *TSTART (2) > 

5 IF  (TSTOPP.GT.'IOOOOO.)  60  TO  9 
GO  TO  (6*7*8) • NSYSOP 

6 TSTOPP«TSTOP 
60  TO  9 

7 TST0PP»TST0PI2) 

60  TO  9 

a TST0PP*AMAX1 (TSTOPIl ) .TST0P(2) ) 

9  CONTINUE 
XSAMEsXSAMEP 

IF  (NSYS0P.EQ.2)  XSAME=XSAHEW 
CALL  TIHEIN  (I*0T*N) 

IF  (NTIN.LT.l)  60  TO  45 
IF  (AREA.LE.O.)  GO  TO  15 
GO  TO  <10*11*12)*  NSYSOP 

10  FLTWl=FLTWTH(l,l) 

FLTW2*FLTV(TH{2*1) 

GO  TO  13 

11  FLTW1*FLTWTH(1*2) 

FLTW23FLTWTH(2,2) 

60  TO  13 

12  FLTWlsAMAXl (FLTWTH (1*1) .FLTWTH (1,2)) 

0FMAXa-FLTW2>RLETH (N) 

13  0FMIN=-FLTW1-RLETH(N) 

0FMAXsFLTW2*RLETH(N) 

0FMIN2»-FLTWTH<1*2)-RLETH(N) 

0FMAX2»-FLTWTH(2*2)*RLETH(N) 

IF  (SYHETRY.GT.O#)  go  to  14 
OFMAXaAMlNl(TAR6ETY^COR»»OTH,OFMAX) 

OFMlNaAMAXl (TARGETY-C0RWDTH*0FMIN) 

14  CONTINUE 
00FS»RLETH (N) /NDOFS 
NOF« (OFMIN-DOFS) /OOFS 

0FMIN«00FS*N0F  t 

NOF> (OFMAX«OOFS) /OOFS 

OFMAXaOOFS*NOF 

15  CONTINUE 

DO  16  NNT«1*NTER 

IF  (TERFR(I}«E0.TFRAC(NNT)«A.HRA0(N)«EQ«ANTH(NNT))  60  TO  17 

16  CONTINUE 
60  TO  45 

17  NALFaNNT 
NNTaNTERA(NALF) 

DO  39  NNTlal*NNT 

ELMIN (1*1) aTERRANE (NNTl * 1 *NALF) ♦CLUTTER 
ELMIN(I*2)>TERRANE(NNTl*2*NALF)*CLUTTER 
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AVSH(1)«AVSH(2)«AVSH(3)*0 
NOFSET ( 1 ) »NOrSET  < 2 ) aNOFSET ( 3 » *0 
OFSET=OFMIN 

IF  < AVERAGE. LE.O.)  OFSETa-SITRAO ( I) 

LAB*10HSAM 

IF  IDEBU61.6T.0.)  PRINT  A6f  LABvFLOAT (N) «FLOAT ( I ) t OFHIN.OFMAX.ELMI 
IN< I *1 ) *ELMIN I 1*2) .float INSTSOP) .TSTARTP^TSTOPP 

18  NSHO*0 
ISHOTaO 
NSR«0 

NSHT ( 1 ) >NSHT ( 2 > >NSHT ( 3 ) >0 

TLW-TSTARTP 

TSaTSTARTP 

T0L«TS 

NPH«1 

IF  (NSYS0P.E0.2)  NPH«3 
NSF»0 

00  36  NT«UNTIN 

IF  (AREA.GT.O.)  GO  TO  21 

QFMIN«THIN(NT*1) 

0FMAXaTHIN(NT*2) 

00FS>IFIX ( ARSIT ( I ) /NOEG) 

OOFSaAHINl (10..00FS) 

00FS»AHAX1 (l.tOOFS) 

IF  (0FHAX>0FMIN.GE.359.9)  OFMAX*OFMAX-.5*OOFS 
0FSET0»0FMIN 

IF  <AVERAOE.LE.O.)  OFSETOaXSITS 

19  XS ITE ( 1) -S ITRAO ( I ) *COS ( OFSETD/OPR ) 

OFSETa-SITRAO ( I ) •S IN ( OFSETO/OPR ) 

NSHT I 1 ) aNSHT ( 2 ) aNSHT < 3 ) a 0 
TLWaTSTARTP 

TSaTSTARTP 

TOL-TS 

NPHal 

NSFaO 

NSHO-0 

IF  (SYHETRY.GT.O.)  60  TO  20 
TEHPOa AMOO (OFSET0.720 • • 360 • ) 

TEMPlaAMOO (TAR6ETY-.5*CORWOTM*720 . *360 . ) 

T£HP2sAM00  CTARGETY^ .S*CORWOTH*720 . *360 . ) 

IF  (TEMP2.GT.TEHP1.A. (TEMPO.lt. TEMPI. 0. TEHPO.GT. TEMP2) ) 60  TO  34 
IF  CTEMP2.LT. TEMPI. A. (TEMPO. LT. TEMPI. A. tempo. 6T.TEMP2))  GO  TO  34 

20  CONTINUE 

21  CONTINUE 
TOaTIN(NT«l) 

0TISaIFlX(.25»TISM(N») 

OTlSaAMAXKOTIS.l.) 

OTalFIXIAMINl I (TIN (NT*2) -TO) /N0TS*.25*T INTER (N) • .25*TINIT (N) ) ) 
OTaAMAXl (0T*1.) 

TOaTO-OT 

TlNTaTlNIT(N) 

22  IF  (NSF.EQ.O)  TOaTO^OT 
IF  (NSF.GT.O)  T0»T0*0TIS 
IF  (TO.GE.TS)  go  to  23 
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NNMT*<TS-TO>/OT^l 

TO«TO*NNNT*OT 

23  CONTINUE 

IF  (FASTRUN.LE.l.)  60  TO  25 
IF  <ISHOT.NE.O>  GO  TO  24 
NSR>NSR-»1 

IF  (NSR.LE.NSOT)  60  TO  25 
IF  (NSF.EQ.O)  TOaTO*NSOT*OT 
IF  (NSF.6T.0)  TO»TO*NSOT*OTID 
60  TO  25 

24  IF  <NSR.6T.NSDT. A. NSF.EQ.O)  TOaTO-NSOT*DT-OT 

IF  (NSR.6T.NS0T.A.NSF.6T.0)  T0»T0-NS0T*0TIS-0TIS 
NSRaO 

25  CONTINUE 

IF  (TO.6T.TIN(NT.2>^.5*0T)  GO  TO  29 
IF  (TS.6E.TIN(NT«2))  60  TO  34 

IF  (0FSET.LT.0FMIN2.0.0FSET.GT.0FMAX2)  60  TO  26 
GO  TO  (26«27«28)*  NSTSOP 

26  Lsl 

GO  TO  31 

27  L»2 

GO  TO  31 
2a  L»1 

IF  (NPH.EQ.4)  GO  TO  31 
IF  (T0.LE.TSTART(2)«TInIT(N))  GO  TO  31 

IF  <T0.6T.TIN<NT»2».O.TO.6T.TST0P(2> .0.TS,GE.TIN1NT,2> .O.TS.GE.TST 
10P(2))  60  TO  29 
L«2 

IF  (NPH.EQ.3)  GO  TO  31 

NPH«2 

L»1 

CALL  INCOV  <I*T0*0FSET,ISH0T»TFIRE*TINT,TAQ*TSTART(2>f2) 

IF  (NSR.GT.NSOT.A.(NSOT>1)*OT*TO,6T.TSTOP(2).A.ISHOT.LE.O.A,FASTRU 
IN.GT.l.)  NSRsO 
IF  (ISHOT.LE.O)  go  to  31 
LAB-IOHGOOD  SHWEA 

IF  (OEBUSl.GT.O.)  PRINT  47*  LAB, FLOAT ( I) *FL0AT (NPH» .FLOAT (NSF) ,TFI 
1RE*T0*TS*T0L*TLW 
IF  (TAQ.LT.TSTART(2))  GO  TO  31 
IF  (NSR.6T.NSDT)  GO  TO  22 
TLW-TO 
NPH«3 
L*2 

IF  (IR(N).NE.O. A. NSF.EQ.O. A. TFIRE.GE.TS)  60  TO  32 
IF  INSF.NE.O)  TS»TS-TISH(N)*TINIT(N) 

IF  (IR(N).NE.O)  GO  TO  31 
TSaAHAXl (TS*T0L^TINIT(N) ) 

60  TO  31 

29  IF  tNSYSOP.NE.3)  GO  TO  34 

IF  (NPH.EQ.l)  GO  TO  34 

IF  (NPH.EQ.2)  GO  TO  30 

IF  (NPH.E0.4)  GO  TO  34 

TO-TLW 

IF  (NSHT(3) .6T.0.A.NSF.6T.0)  TS=TS-TISH(N) ♦TINIT (N) 
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IF  (IR(N).EQ.O)  TS«AMAXl<TS«TOL«T1NIT(N)) 

LAB«10HCHAN6E  1 

IF  (OEBUttl.GT.O.)  PRINT  <»7t  LABtTS»TO*AVSH  (3) 

30  L*1 
NPH«4 

31  CALL  INCOV  4 1 •TO»OFSET» ISHOT»TFIRE*TINT«TAQf TStL) 

32  IF  (ISHOT.LE.O)  GO  TO  22 
IF  (NSR.GT.NSOT)  GO  TO  22 
TOL-TO 

IF  (AREA.GT.O.)  OFSETOsOFSET 
LAB-IOHGOOO  SHOT 

IF  tOEBUGl.GT.O.)  PRINT  ^7*  LABtOFSETOt TS«TAQ« AZTA*ELTA*RNGTA« ASPT 
lA*TFIRE*AZTF*ELTF.RNGTF*ASPTF.TOf AZT*ELT*RNGT*ASPTtNPH 
IF  <NPH»NE.3»A.NSYSOP.NE.2)  AVSH (2) sAVSH (2> • 

IF  (NPH.NE»3.A.NSYS0P.NE»2)  NSHT(2)=1 

IF  (NPH.NE.3.A.NSYS0P.NE.2.A.TAQ.GT.TEGRESS)  AVSM(2)*AVSH(2) ♦FEGRE 
1SS«1 . 

IF  (NPH»NC.3.A.NSYS0P.NE.2.A.T0.LE.TSTART(2)1  AVSH ( 1) sAVSHi 1 ) ♦! . 

IF  (NPH»NC.3«A.NSYS0P.NE*2.A.T0.LC.TSTART(2> I NSHT<1)«1 
IF  (NPH.NE.3.A.NSYS0P.NE.2.A.T0.LE.TSTART (2) .A.TAQ.GT.TEGRESS)  AVS 
1H41)*AVSH(11^FEGRESS-1. 

IF  (NSYS0P,NE.1,A.NPH,E0.3,A,TAQ,GT.TSTART(2>)  AVSH(3)sAVSH(3)*1. 
IF  (NSYS0P.NE.l.A.NPH.EQ*3.A.TAQ.GT.TSTART(2n  NSMT(3»»1 
NSF«NSF^1 
TS*TFIRE*TISHIN) 

IF  (NSF.LT.NSS IN))  60  TO  22 
NSF«0 

TS»TFIRE*TINTER|N) 

NSH0«NSH0«1 

IF  (TS.6E.T0«TInIT(N))  GO  TO  33 
TS«TO^TINIT(N) 

33  CONTINUE 

IF  (NSHO.LT.NS(N) ) GO  TO  22 
NSHO*0 

TS«TFIRE*TRELOAO(N) 

IF  (TS.6E.T0«TINIT(N) ) GO  TO  22 
TS-TO*TINIT(N) 

GO  TO  22 

34  CONTINUE 

IF  (AREA.GT.O.)  GO  TO  36 
DO  35  K«l«3 

35  IF  (NSHT(K).GT.O)  NOFSET (K) *NOFSET IK) ♦! 

LAB>10HPOLARLIMS 

IF  (OEBU62.GT.O.)  PRINT  46*  LABtFLOAT INT) tOOfSyOFSETDtXSITE ( I ) .OFS 
lET*AVSM.FLOAT(NOFSET)fTSTART 
IF  (AVERAGE. LE.O.)  60  TO  36 
OFSETO*OFSETD^DOFS 
IF  (OFSETO.LE.OFMAX)  80  TO  19 

36  CONTINUE 

IF  (AREA.LE.O.)  GO  TO  38 
00  37  Kal,3 

37  IF  (NSHT(K).6T.O)  NOFSET (K) ■NOFSET (K> ♦! 

IF  (AVERAGE. LE.O*)  60  TO  38 
OFSETsOFSET^OOFS 
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IF  (OFSET.LE.OFMAX)  60  TO  16 

38  CONTINUE 

DO  39  K«l,3 

IF  (NOFSETIKI.GT.O)  AVSH (K) sAVSH (K) /NOFSFT (K ) 
SITW0TH(I«K>=N0FSET(K)*00FS*TERRANE(NNT1«3*NALF)>SITW0TH(I,K) 
AVSMOT ( I tK ) sAVSHOT ( I vK ) ♦TERRANE (NNT 1 » 3»NALF ) •AVSH (K ) 

LABsIOHAVSHOTS 

IF  <OEBU61.GT.O.)  PRINT  46*  LAB*FLOAT ( I ) tFLOAT (NOFSET (K) ) , AVSHOT ( I 
1»K)*SITW0TH(I»K) tAVSH(K)*NOFSET(K)»OOFS 

39  CONTINUE 
XSITE(I>«XSITS 
00  40  Kalf3 

40  IF  (AVERAGE.LE.O.)  SITWOTHI I*K)*CORWOTH 

IF  (I.EQ.NTOTS.O.AREA.LE.O..O.XSAME.EO.O.)  60  TO  45 
K»1 

IF  (NSYS0P.EQ.2)  K=2 
XOONEaXSAHE^XSTART-RLETHCN) 

IF  <XSITE(I).GT.XOONE)  60  TO  45 

XTSTOPxTRP (TSTOPtT  C 1 »K ) tXYZ ( 1 • 1 1 K ) »NXYZT (K) ) ♦XSTART 
XTSTART»TRP<TSTART*Tn*K)  fXYZ(l*l*K)  .NXYZT(K)  )*XSTART 
XTE6R=TRP (TEGRESStT ( 1 f K ) *XYZ ( 1 • 1 «K ) tNXYZT (K » ) ♦XSTART 
IF  (ABS(XTEGR-XS1TE(I) ) •LT.RLETH(N) .A.FEGRESS.NE.I.)  GO  TO  45 
D31aXTEGR-XSITE(I) 

IF  (NSYS0P.EQ.2)  60  TO  41 

XTREL*TRPITSTART<2) ♦T<ltK) *XYZ( 1.1 *K) *NXYZT(K) » ♦XSTART 
IF  (ABS(XTREL-XSITE(I) ) .LT.RLETH(N) ) 60  TO  45 
D4I*XTREL«XSITE(I) 

41  CONTINUE 

IF  ( ABS  CXTSTART-XSITE ( I ) ) .LT .RLETH (N) .0. ABS ( XTSTOP-XS ITE ( I ) ) .LT.RL 
lETHINl)  60  TO  45 
OlIaXTSTART-XSITEd) 

02I»XTST0P-XSITE(n 

J1»M1 

00  44  JsJlfNTOTS 

IF  (NTYPECII .NE.NTYPEIJ) .O.TEPFR ( I ) .NE.TERFR ( J) .O.XSITE( J) .GT.XnON 
lE.O.NOONE(J) .GT.O)  GO  TO  44 

IF  <ABS(XTSTART-XSITE< J> ) .LT.PLETH<N> .O.ABS(XTSTOP-XSITE( J) ) .LT.RL 
1ETH(N)>  60  TO  44 
OlJaXTSTART-XSITE { J) 

02 JaXTSTOP-XSITE  < J) 

03J>XTEGR-XSITE(J) 

IF  (FEGRESS.NE.1..A.S6N(D3I).NE.S6N(D3J))  60  TO  44 

IF  (NSYS0P.EQ.2)  60  TO  42 

04J«XTREL-XSITE(J» 

IF  (SGN<D4I) .NE.S6N(04J) ) GO  TO  44 

42  CONTINUE 

IF  (SON(DII) .NE.SGNIDlJ) .0.S6N(02I) .NE.SGN(02J) ) GO  TO  44 
00  43  Kal,3 

SITWOTH(J,K)aSITWOTM(I,K) 

AVSH0T<J.K)aAVSH0T(I.K) 

43  CONTINUE 
NOONEf J>sl 

44  CONTINUE 

45  CONTINUE 
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TSTARTPaSVTl 

TST0PP»SVT2 

RETURN 


46  FORMAT  llXvAlOtlZFlO.A) 

47  FORMAT  ( 1X*A10« 17F7.1 » I3) 
END 


nn  on  nn 


SUBROUTINE  PROBS 


1 


THIS  ROUTINE  CALCULATES  AND  PRINTS  PROBABILITY  OF  SURVIVAL  FOR 
EACH  XSTART 


COMMON  /TRAJ/  NXY2T(2) ,XYZ ( 1000*3*2) *T ( 1000*2) »XSTART*TST ART (2) tTS 
1T0P(2) *T IN (10*2) yNTINtAECM (20*2) *NXECM,THIN(1 0*2) *0XSTART*NXSTART* 
2XSAMEP*XSAMEH*NSYS0P*TST0PP*TSTARTP*TARGETY*0XYZ0T( 1000*3*2)* IS (2) 
COMMON  /ISAM/  NTYPE ( lOO) *NSITE ( 100 ) *XSITE ( 100 ) *PUP ( 100 ) *ELMIN( 100* 
12) *NTOTS*SIT)iDTH(100*3)*AVSHOT(100*3)*rERFR(100)*SITRAD(100)*ARSIT 
2(100) 

COMMON  /AS AM/  HRAO ( 1 0 ) *RTRK ( 1 0 ) * ELM AX (10)*TINIT(10)*TINTER(10)*NS( 
110>*NSS(10)*TRElOAO(10)*AVVEL(10)*ASPMIN(10)*ASPMAX(10)*AZMAX(10)* 
2RN6(20*10)*ELR(20*10)*NRN6(10)*FUS(20*10)*ELF(20*10)*NFUS(10>*IP(1 
30)*TISH(10)*ECME(10)*AlTMIN(10)*ALTMAX(10)*SI6TH(20*4)*SIG(20*4)*N 
4SI6(4)*RLOCK(10)*XM1,SL(20«10)*THISL(20*1C)*NXHISL(10)*RAOTRK(10)*R 
SAOLOCKdO) 

COMMON  /PARM/  DF(IO) *0EauG*FLTHTH(2*2) *RLETH(10) *DEBUG1 *TITLE (6) *P 
1KSS(10*2)*CORWOTH*CLUTTER*TERRANE<10*3*10)*NTERA(10)*NTER*TFRaC(10 
2) *ANTH(10)*AREA,OEBUG2*OFSETO*RELEASE*VELPCN*OLEV(10) *ndlev*averag 
3E*NITLI (6)*NITL2(6)*SYMETRY*TEGRESS*FEGRESS*NITL3(6) *NINrRPR ( 10) *N 
4ASP*FASTRUN 

COMMON  /AVG/  SUMECM ( 10*3) «SUMNECM( 10*3) 

DIMENSION  PECM(10*3>*  PNECH(10*3) 

PRINT  14 

IF  (AREA.LE.O.)  PRINT  17 
IF  (AREA.LT.O.)  PRINT  18 
IF  (AVERAGE.GT.O.)  PRINT  15 
PRINT  16*  XSTART 
00  1 1^1*30 

1 PECM(I)>PNECM(I)sl. 

DO  12  1*1*3 

11*1 

IF  (I.EQ.3)  11*2 
00  12  J*l*NTOTS 
N*NTYPE(J) 

IF  (DF(N).LE.O.)  go  to  12 
ENC-1. 

IF  (AVERAGE.LE.O.)  go  to  6 
IF  (SYMETRY.LE.O.)  go  to  2 
ENC*S I T WOTH ( J • I ) /COR MOTH 
GO  TO  6 

2 IF  (AREA.LE.O.)  GO  TO  4 

3 ENC*SITH0TH ( J* I )/CORWOTH- ( AMINl ( 1 . * .5*SITW0TH( J* I ) /CORWOTH) ) *•2 
GO  TO  5 

4 IF  (SITWOTH(J*I) .GE.360..0.SITMDTH(J*I)/CORWOTH.GC.2.)  GO  TO  6 
IF  (CORWOTH*. S»SITMOTH(J*I) .LE. 360.)  GO  TO  3 
Al*C0RW0TH*.S*SlTW0TH(J*I)-360. 

A2*C0RW0TH-.5*SITW0TH(J*I) 

PI*  (COR*(OTH*SITWOTH(J»l) -360.) /CORWOTH 

ENC* (SITWOTH ( J* I ) ♦ ( A 1 *A2)* (Pl-S ITWOTH ( J, I ) /CORWOTH) ) /CORWOTH 

5 ENC*AHIN1 (l.vENC) 


42 


n o 


r 


6  EI*NSITEIJ>*ENC 

Sl«El*AVSHOT(J.n*PUPIJ)*Of  (Nl 

ECMEFF*! • 

IF  <NXECM.LC*0)  go  to  8 

??I,Sm<5u6E.XECH,L.l..*.«MTE.J..LE.XECH,t.2..  ECHEEr-ECHE.N. 

7 CONTINUE 

8 CONTINUE 

00  11  L*1*N0LEV  « o 

rei.<i.-P«ss<N.ii)>"<si‘i)>.Ev(L)) 

9 “tKl‘.ipKSSIN.nfECHEFr.SI-OLEV.L> 

10  CONTINUE  e 

PECNIL*1>*^ECM<L*  I1*I^SIE 

PNECMIL*I)*PNECM<L*I1*PSI 

11  CONTINUE 

12  CONTINUE 

.0LEV.U.<PNECH,L.n.....3...PEC«.L.....>..3..PFCH.L...- 

lPNEci<L.3>.L*l.N0LEV) 

SwIJ.  ,NTTPE(l..KSnE.U.SI«.0,I)..SIT«TmI.J...VSHOT<I.J>.J 

l*If3)*I»l*NT0TSl 
00  13  1*1*3 

SShKh'JlI  U^KeCH  1L. I ) -PECH (L . I ) '>!« 

. I ) -SUHNECH  a . n *pnecm  il.  i . /nxsiapt 

13  CONTINUE 
RETURN 


U JSSS"t  !ls;'«s”}r;iU.OED  over  .LL  OPPSETS  .PPL.CA8LE. 

}‘t  !!6H’'II*or.L%25;viL. 

inrti 

v/o  ech. 

20  FORMAT  (8F1S.31  a*.AHY  OR  R*3(IX*9HLETH  WDTH*2X*8HA 

COHPLETE.-.X..9H-EAPOR  . 

2REL  TO  TAR)) 

22  FORMAT  (I5*0F1O«3) 

ENO 
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SUBROUTINE  INPUTS 


C 

C 

C THIS  ROUTINE  READS  ALL  INPUT  DATA  AND  TERMINATES  THE  PROGRAM 

C 

C 

COMMON  /TRAJ/  NXYZT(2) «XYZ ( 1000*3«2) >T < 1000*2) tXSTARTtTSTART <2) fTS 
1TOP(2>*TIN(10*2)*NTIN*XECH(20«2)*NXECH*THIN(10«2)»OXSTART«NXSTART* 
2XSAMEP*XSAMEVltNSYSOP*TSTOPP*TSTARTP«TARGETY*OXYZOT(1000«3*2)«IS(2) 
COMMON  /ISAM/  NTYPE ( 100) •NSITE ( 100) *XSITE ( 100) •PUP( 100) *ELMIN( 100* 
12) *NTOTS«SITWOTH 1100*3) yAVSHOTi 100*3) *TERFR(100) *SITRAO<100) *ARSIT 
21100) 

COMMON  /ASAH/  HRAO (1 0 ) *RTRK ( 1 0 ) • ELMAX (10)*TINIT(10)*  T INTER ( 1 0 ) * NS ( 
110)*NSS(10)*TRELOAO(10)*AVVEL(10)*ASPMIN(10)*ASPMAX(10)*AZMAX(10)* 
2RN6(20*10)*ELR(20*10)*NRNG(10) *FUS (20* 10) *ELF ( 20* 10 ) *NFUS ( 10) * IR < 1 
30)*TISH(10)*ECME(10)*AlTMIN(10) *ALTMAX(10>*SIGTH(20*4)«SIG(20*4) *N 
4SIG(4)*RLOCK(10)*XHISL(20*10)*TMISL(20*10)*NXMISL(10) *RAOTRK(10)*R 
SAOLOCKdO) 

COMMON  /PARM/  OF(IO) *0EBUG*FLTWTH(2*2) *RLETH(10) *0EBUG1*TITLE(6) *P 
1KSS(10*2)*CORWOTH*CLUTTER*TERRANE(10*3*10)*NTERA(10)*NTER«TFRAC(10 
2) *ANTH(10)*AREA*OEBU62*OFSETO*RELEASE*VELPEN,OLEV(10)*NOLEV*AVERAG 
3E*NITL1(6)*N1TL2(6)*SYMETRY*TEGRESS*FEGRESS*NITL3(6)«NINTRPR(10)*N 
4ASP*FASTRUN 
DIMENSION  NTEP(6) 

INTEGER  TITLE  * 

DATA  OEBUG/0./*OEBU61/0./*OF/10*0./*NTOTS/0/*NXYZT/2*0/*XSTART/0./ 
1*TSTART/2«0./*TSTOP/2«0./*NXECM/0/*CLUTT£R/.2S/*NTER/0/*AREA/1./,0 
2EBUG2/0./ 

DATA  DXSTART/2.5/*NXSTART/1/ 

DATA  DLEV/1,*9*0./*NDLEV/1/*NSIG/4*0/*PKSS/20*0./*TSTARTP*TSTOPP/2 
1*-100000./ 

DATA  AYERA6E/1./ 

DATA  SYMETRY/1./ 

DATA  FASTRUN/0./ 

DATA  NITLI/12*10H  /.NITL3/6*lOH  / 

DATA  TARGETY/0./ 

DATA  TEGRESS/100000./*FEGRESS/1./ 

DATA  RAOTRK/10*0«/*RADlOCK/10*0./*NINTRPR/10*1/ 

IER>0 

DATA  TITLE/6*10H  / 

DATA  NITL2/6*10H  / 

PRINT  60 

1 READ  61*  N*X*NTEP 
PRINT  62*  N*X*NTEP 

IF  (N.NE.7HENDCASE)  GO  TO  2 
IF  (lER.EO.O)  return 
STOP 

2 IF  (N.EQ.6HENDJ0B)  STOP 
IF  (N.NE.SHTITLE)  00  TO  4 
DO  3 J»l*6 

3 TITLE(J)=NTEP(J) 

GO  TO  1 

4 IF  (N.NE.4HXYZT)  GO  TO  8 
DO  5 J*l*6 
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5 NXTL2<J»*NTEP<J) 

XSAMCPaX 

PRINT  63 
X>NXYZTIl)sO 

6 READ  64*  <XVZ(K«l*J«l)f J31,3)9T(K«1«1) 

IF  (EOFdl.NE.O)  GO  TO  1 

PRINT  64t  (XYZ(K*l»J»l)»J«l>3)»T(K«ltl) 

NXYZT<l>*K«K^l 
IF  (K.EQ.l)  GO  TO  6 
DO  7 J*l*3 

7 OXYZOT<K#J»l)«<XYZfK»J,I)-XVZIK-l*J»l)»/(T(K*l)-T(K-lfl>) 

GO  TO  6 

8 IF  (N.NE.6HTSTART)  GO  TO  9 
TSTART»X 

GO  TO  1 

9 IF  (N.NE.6HXSTART)  GO  TO  10 
XSTART»X 

GO  TO  1 

10  IF  (N.NE.SHTSTOP)  GO  TO  11 
TSTOPsX 

GO  TO  1 

11  IF  (N,NE.4HSITE)  GO  TO  14 
DO  12  J»l*6 

12  NITLl f J><NTEP(J) 

PRINT  65 
NTOTSaO 

13  READ  66*  NTYPE<NTOTS^l ) fNSITE (NTOTS^l ) tXSITE INTOTS*! ) *SITRAO (NTOTS 
l*l)*PUPINTOTS*l)*TERFR(NTOTS*l) 

IF  (EOFdl.NE.O)  GO  TO  1 

PR INT  66*  NTYPE (NTOTS* 1 ) «NSITE (NTOTS* 1 » * XSITE  <NTOTS* 1 » *5 ITRAO (NTOT 
1S>1»*PUP(NTOTS>1I*TERFR(NTOTS*1) 

NTOTSsNTOTS*! 

GO  TO  13 

14  IF  (N.NE.3HSAH)  GO  TO  15 
J«X 

READ  67*  NSfJ) *NSS ( J) • IP ( J) «HRAO(J) *RTRK ( J) »RLOCK ( J) tELMAX < J) , ALTH 
1IN(JI 
PRINT  68 

PRINT  67*  NS(J) «NSS(J) vIRlJ) »HRAO(J) «RTRK ( *RLOCK (J) *ELHAX ( J) * ALT 
IMIN(J) 

READ  64*  TINIT(J>*TISH(J)tTINTER(J)*TRELOAO(J>*ASPMIN(J)«ASPMAX(J) 
1*AZMAX(3)«ECME(J) 

PRINT  69 

PRINT  64.  TINIT(J) .TISH( J) .TINTERI J) .TRELOAOl J) .ASPMINI J) .ASPMAX ( J 
1).AZMAX(J).ECME(J) 

PRINT  70 

READ  64.  RAOTRK(J) .PADLOCK Ul 
PRINT  64.  RADTRKf J).RAOLOCK(J) 

GO  TO  1 

15  IF  fN.NE.3HRN6>  GO  TO  17 
K«0 

J*X 

RLETH<J)«>1.E*10 

ALTNAX(J)»-1.C*10 
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PRINT  71 

16  READ  66*  ELRtK«l*J)*RN6(K*l»J) 

IF  (eof(I).ne.o)  go  to  1 
PRINT  64»  ELR(K«lf J)fRNG(K«lf J) 

NRNG(J)sK«K«l 

RLETH(J)«ANAX1 fRLETHIJ) •RN6<K> J) •COSlELR fK* J)/57.2958) ) 
ALTHAX(J)«AMAX1  (ALTMAXI  J>  •RNG(K«J)«SIN(ELR(K,JI/S7.2<158) ) 
GO  TO  16 

17  IF  (N.NE.AHFUSE)  GO  TO  19 
K»0 

J>X 

PRINT  72 

18  read  66*  ELF(K4l»J>tFUS(K»ltJ) 

IF  (EOF(1),NE.O)  GO  TO  1 
PRINT  64»  ELF(K«ltJ)«FUS(K«l«J) 

NFUS(J)sKsK^l 

GO  TO  18 

19  CONTINUE 

IF  (N.NE.2H0F)  GO  TO  22 
IF  (X.LE.O.)  GO  TO  21 
DO  20  1«1«10 

20  DF(I)«X 
GO  TO  1 

21  READ  73*  OF 
PRINT  73 t OF 
GO  TO  1 

22  CONTINUE 

IF  (N.NE.SHOEBUG)  GO  TO  23 

OEBUGsX 

GO  TO  1 

23  CONTINUE 

IF  (N.NE.6HOEBUG1)  GO  TO  29 
0EBUG1«X 
GO  TO  1 
29  CONTINUE 

IF  fN.NE.THCORMDTH)  GO  TO  25 

CORWDTHaX 

GO  TO  1 

25  IF  (N.NE.9HPKSS)  GO  TO  26 
J»1 

IF  (X.NE.O.)  J32 

READ  73t  (PKSS(K»J)»Ksi,10) 

PRINT  73*  (PKSS(K»J)«Ksl*10) 

GO  TO  1 

26  IF  (N.NE.9HXECM)  GO  TO  28 
J»0 

NXECMaO 

27  READ  69«  XECM ( J«1 t 1 ) fXECM IJ^l •2) 

IF  (EOF(l).NE.O)  GO  TO  1 

PRINT  69*  XECM(J«l«l)tXECM(J»lf2) 

NXECMaJaJ^l 
GO  TO  27 

28  CONTINUE 

IF  (N.NE.7HCLUTTER)  60  TO  29 
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CLOTTER-X 
GO  TO  1 

29  IF  <N,NE«7HTERRAIN)  GO  TO  32 
NTER-0 

30  READ  64.  TFRACCNTER*! ) , ANTHINTER+l ) 
IF  (EOF(l).NE.O)  60  TO  1 

PRINT  74 

print  64.  TFRACINTER*!) »ANTH(NTER>1) 
NTER»NTER*1 


J«0 

NTERA INTER) sO 
PRINT  7S 

31  READ  64.  TERRANECJ>l»l.NTER)*TERRANEU*1.2.NTER)*TERRANE<J*lf3*NTE 

IR) 

IF  (EOFID.NE.O)  60  TO  30 

PRINT  64.  TERRANEl J»1 .1 .NTER) .TERRANE ( J*l .Z.NTER) .TERRANE IJ* 1 .3.NT 
lER) 

MTERA<NTER)*J*J>1 
60  TO  31 

32  CONTINUE 

IF  (N.NE.4HAREA)  60  TO  33 


AREA«X 
60  TO  1 

33  CONTINUE 

IF  (N.NE.6H0EBUC2>  60  TO  34 

DE8U623X 

GO  TO  1 

34  CONTINUE 

IF  IN.NE.TMTSTARTW)  60  TO  35 

TSTARTI2)=X 

60  TO  I 

35  CONTINUE 

If  (N.NE.6HTST0P«)  GO  TO  36 

TST0P<2)»X 

60  TO  1 

36  CONTINUE 

IF  IN.NE.7H0XSTART)  GO  TO  37 

OXSTART>X 

GO  TO  1 

37  IF  (N.NE.7MNXSTART)  60  TO  38 

nxstart«x 

GO  TO  1 

38  CONTINUE 

IF  (N.NE.7HAVERA6E)  60  TO  39 

AVERAGEsX 

60  TO  I 

39  CONTINUE 

IF  (N.NE.4H0LEV)  GO  TO  41 
NOLEVsJ»0 

40  READ  64.  0L£V(J*1) 

IF  lEOFfD.NE.O)  60  TO  1 
PRINT  64.  OLEVIJ^l) 
N0LEV*J*J*1 
GO  TO  40 
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41  CONTINUE 

IF  tN*NE.7HSYMETRY)  GO  TO  42 

symetry*x 

60  TO  1 

42  CONTINUE 

IF  (N.NE.6HE6RESS)  60  TO  43 
READ  64«  TEGRESS*FEGRESS 
PRINT  76 

PRINT  64.  TEGRESS.FEGRESS 
GO  TO  1 

43  CONTINUE 

IF  (N.NE.SHXYZTW)  60  TO  48 
00  44  J«1.6 

44  N1TL3<J»*NTEP{J) 

XSAMEW>X 

PRINT  63 
K«0 

45  READ  64.  (XYZ (K.l • J.2> * J=1 .3) *T (K«l *2) 

IF  (EOF(IUNE.O)  GO  TO  47 

PRINT  64.  (XYZ<K«l.J«2).Jsl.3).T(K>1.2) 

KaK.l 

IF  <K.EQ.l>  GO  TO  45 
00  46  Jal.3 

46  0XYZ0T(K.J,2»*(XYZ<K.J,2)-XYZ<K-1.J,2))/(T<K.2)-T(K-1.2) ) 
GO  TO  45 

47  NXYZT<2>«K 
GO  TO  1 

48  IF  (N.NE.9NSIGNATURE)  GO  TO  51 
J»X 

NS«0 

PRINT  77 

49  READ  64.  SIGTH(nS«1 . J) «SI6(NS*1 . J) 

IF  (EOF(I).NE.O)  go  TO  50 

PRINT  64.  SI6TH<NS«1.J) .SIG(NS.1«J) 

NSaNS«l 
GO  TO  49 

50  NSIGUIaNS 
GO  TO  1 

51  IF  (N.NE.TSTOPW)  GO  TO  52 
TST0P(2)aX 

GO  TO  1 

52  CONTINUE 

IF  (N.NE.6HMISLXT)  GO  TO  54 

PRINT  78 

JaX 

NXMISL(J>aN*0 

53  READ  64.  XNISL fK*l . J> .THISL (K.l . J) 

IF  (EOF(IUNE.O)  60  TO  1 

PRINT  64.  XMISL(K«1»J>,TMISL(K.1,J) 

NXMISL(J)«K«K«I 
60  TO  53 

54  IF  (N.NE.6HTST0PP>  GO  TO  55 
TSTOPP-X 

GO  TO  1 
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55  IF  (N.NE.7HTSTARTP)  GO  TO  56 
TSTARTPsX 

GO  TO  1 

56  CONTINUE 

IF  (N.NE.5HN0SIG)  GO  TO  57 
READ  79t  NINTRPR 
PRINT  79*  NINTRPR 
GO  TO  1 

57  CONTINUE 

IF  {N*NE.7HFASTRUN>  GO  TO  58 

FASTRUN«X 

GO  TO  1 

58  CONTINUE 

IF  (N.NE.TARGETY)  GO  TO  59 

TARGETY»X 

GO  TO  1 

59  CONTINUE 
PRINT  80 
lERal 

GO  TO  1 
C 


60  FORMAT 

61  FORMAT 

62  FORMAT 

63  FORMAT 

64  FORMAT 

65  FORMAT 
ITERFR) 

66  FORMAT 

67  FORMAT 

68  FORMAT 


(IHI) 

(A10«F10.0*6A10) 

(1X*A10*F10.4*10X*6A10) 

(9X*1HX*9X«1HY»9X*1HZ*9X*1HT) 

(8E10.3) 

(SX*5HNTYPE*5X*5HNSITE*5X*5HXSITE*5X*5HYSITE*7X.3HPUP*5X,5H 


(2I10*6F10.4) 

(3I10t5£l0*3) 

( 8X • 2HNS  * 7X  * 3HNSS  * SX  * 2H I R « 6X  » 4HHR  AO • 6X • 4HRTRK  « 5X  * 5HRL0CK  « 5X 
1»5HELMAX*4X*6HALTMIN) 

69  FORMAT  C5X.5HTINIT*6X*4HTISM*4X*6HTINTER*3X«7HTREL0A0*4X*6HASPMIN* 
14X  *6HASPMAX*5X  vSHAZMAX ,6X  «4HECME ) 


70  FORMAT  (4X*6HRA0TRK«3X*7HRA0L0CK) 

71  FORMAT  (7X*3HELR«7X»3HRNG) 

72  FORMAT  {7X*3HELF*7X*3HFUS) 

73  FORMAT  (10F8.2) 

74  FORMAT  (SX*SHTFRAC»6X*4HANTH) 

75  FORMAT  I4X«6HAV  ELF«4X,6HAV  ELR«6X*4HPR0B) 

76  FORMAT  <3X.7HTEGRESS.3x.7HFE6RESS> 

77  FORMAT  (5X«5HSIGTH*7X»3HSIG) 

78  FORMAT  (5X«5HXMISL»5X*5HTMISL) 

79  FORMAT  (1018) 

80  FORMAT  (46H  PREVIOUS  CARD  NOT  IDENTIFIED*  JOB  TERMINATED  ) 
END 
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SUBROUTINE  IMCOV  <NSAM,TI .OFSET* lSHOT»TFIRE*TlNT.TAQf TS.Ll ) 


THIS  ROUTINE  DETERMINES  IF  A SAM  CAN  SHOOT  A MISSILE  TO  INTERCEPT 
AT  TIME  TI 


TI  TIME  AT  intercept 

NSAM  INDEX  NUMBER  OF  SAM  TO  BE  CONSIDERED 

OFSET  OFFSET  DISTANCE 

TFIRE  TIME  SAM  WOULD  HAVE  HAD  TO  FIRE  FOR  INTERCEPT  AT  TI 

ISHOT  ■!  CAN  GET  OFF  A SHOT  WITH  INTERCEPT  AT  TI 
COMMON  /TRAJ/  NXYZT (Z) ,XYZ ( 1000*3*2) *T (1000*2) *XSTART*TSTART (2) *TS 
lTOP(2)*TlN(10*2)*NriN*xECH(20*2)«NXECM*THlN(10*2)*OXSTART*NX5TART* 
2XSAMEP*XSAMEW*NSYSOP*TSTOPP*TSTARTP*TAR6ETY*OXYZOT(1000*3*2)*IS(2) 
COMMON  /ISAM/  NTYPE(IOO) *NSITE < lOO) *XSITE ( 100) *PUP(100) *ELMIN(100* 
12) *NT0TS*SITW0TH( 100*3) *AVSH0T( 100*3) *TERFR( 100) *5 ITRAD( 100) *ARSIT 
2(100) 

COMMON  /ASAM/  HRAO(IO) ,RTRK(10) *ELMAX(10) *T INIT ( 10 ) *TINTER ( 10) *NS ( 
110)*NSS(10>*TRELOAO(10)*AVVEL(10) *ASPMIN(10) *ASPMAX(10) •AZMAX(IO) * 
2RNG(20*10)*ELR(20*10)*NRN6(10)*FUS(20*10)*ELF(20*10)*NFUS(10) *IR(1 
30) *TISH(10)*ECME(10) *AlTMIN(10) *ALTMAX(10) *SIGTH(20*4) *SIG(20*4) *N 
4SIG(4) *RLOCK(10) *XMISL(20*10) *TMISL(20*10) *NxMISL(10) *RA0TRK (10) «R 
SADLOCKdO) 

COMMON  /P ARM/  OF ( 1 0 ) *OEBU6 * FLT WTH ( 2 * 2 ) « RLETH (10)*  DEBUG 1 « T I TLE ( 6 ) * P 
1KSS(10*2) »CORWOTH*CLUTTER*TERRANE(10*3*10) *NTERA(10) *NTER*TFRaC(10 
2) *ANTH(10)*AREA,OEBUG2,OFSETD*RELEASE*VELPEN,DLEV(10)*NOLEV*AVERAG 
3E*NITL1(6)*NITL2(6)*SYmETRY*TEGRESS*FEGRESS*N1TL3(6)*NINTRPR(10)*N 
4ASP*FASTRUN 

COMMON  /o/  azt*elt*rngt*aspt*aztf*eltf*rngtf*asptf*azta*elta*rngta 

1*ASPTA 
K«Ll 
ISHOTaO 

NSTaNTYPE(NSAM) 

IF  (OF(NST).LE.O.)  RETURN 
C FIND  AZ*EL*ASP*RNG  OF  PENETRATOR  WRT  SITE  AT  TI 
NASP-1 

CALL  COORD  (NST*XSITE(NSAM) *OFSET*Tl *AZT*ELT.RN6T*ASPT*K) 

IF  (RLETM(NST).LT.RNGT*C0S(ELT/57.2985))  RETURN 
ALTaRNGT*SIN(ELT/57,29sa) 

IF  (ALT.LT.ALTMIN(NST))  RETURN 
NASP*0 

C CALCULATE  FLYOUT  TIME*FIRING  TIME*  ACQUISITION  TIME 

TFIRE*TI-TRP(RNGT*XMISL(1*NST)*TMISL(1*NST)*NXMISL(NST) ) 

TAO«TF IRE-TINT 
NN>10HTFIRE*TAQ 

IF  (OEBUG.GT.O.)  PRINT  12*  NN*TFlRE*TAO*TI *TS*TSTART 
IF  (TFIRE.lt. TS)  RETURN 

IF  (TFIRE. LE.TSTARTP.O.TAO.LT.TSTARTP)  RETURN 
C CHECK  IF  WITHIN  LETHAL  ENVELOPE 
NN«I0H  ENVEL  CK 

IF  (OEBUG.GT.O.)  PRINT  12*  NN*TI*OFSET*AZT*ELT*RNGT*ASPT*RMIN*RMAX 
l*FLOAT(K) 

IF  (AZT.6T.AZMAX(NST))  RETURN 
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L«1 

IF  (AZT.8T,90.)  L«2 

IF  (IR(NST).LE.0*A«(CLT*LT.ELMIN(NSAH«L) .0.ELT.GT.CLHAX<NST) ) ) PET 
lURN 

RMAX*TRP (ELT  tELR ( 1 tNST) »RNG ( 1 f NST) «NRNG (NST ) ) 

IF  (RNGT.6T*RMAX)  return 

RM IN«TRP ( ELT  t ELF ( 1 • NST ) « FUS ( 1 • NST ) t NFUS ( NST ) ) 

IF  (Rngt.lt»rhin)  return 

C FIND  AZ*EL«RN6>ASP  OF  PENETRATOR  NRT  SITE  AT  TFIRE 
K»L1 

IF  (NSYSOP.Ea.3.A.TFIRE,LT.TSTART(2).A.K.E<3*2)  K»1 
K1>K 

IF  (IRfNSTl.GT.O)  K1*K1«2 

CALL  COORD  <NST,XSITE(NSAM) *OFSET*TFIRE*AZTF,ELTF»RNGTF*ASPTF,K) 

IF  (NINTRPR(NST) .GT.O)  GO  TO  1 

APENaTRP<ASPTF,SlGTH(l,Kn.SIG(l,Kl).NSlG(KU) 
RMTF»SQRT(APEN/RL0CK(NST> ) 

IF(IR<NST).LE.O)RMTF=SQRT(RMTF) 

IF  fNINTRPR(NST).LT,0)  GO  TO  2 
GO  TO  3 

1 RMTF=RAOLOCKINST) 

2 IF  (ASPTF.LT.ASPNIN(NST)*0.ASPTF.6T.ASPMAX(NST))  RETURN 

3 CONTINUE 
NNbIOHFIRE  CK 

IF  (OEBUO.GT.O.)  PRINT  12»  NN»AZTF»ELTF#RN6TF* ASPTF,FLOAT CK) *FP.OAT 
ICKI),RNTF 

C CHECK  IF  WITHIN  FIRING  CONSTRAINTS  AT  TFIRE 
L»1 

IF  (AZTF.GT.90.>  L*2 

IF  (RNGTF.GT.RHTF.O.ELTF.LT.ELHINCNSAMvL) .O.ELTF.GT.ELMAXCNST) ) PE 
ITURN 

C CHECK  IF  OBSERVED  AT  TAQ 
K»L1 

IF  (NSYS0P.EQ.3.A.TAQ.LT.TSTART(2).A.K.EQ.2)  K=1 
K1«K 

IF  (IR(NST) .GT.O)  K1=K1*2 

CALL  COORD  (NST.XSITE CNSAM) *OFSET*TAO. AZTA»ELTA»RNGTA. ASPTA.K ) 

IF  CNINTRPR<NST).GT.O)  GO  TO  4 

APEN«TRPCASPTA,SIGTH<1,K1).SIG<1,K1) *NSIG(K1) ) 

RNTA-SQRT  C APEN/RTRK (NST) ) 

IF (IRCNST) .LE.OIRMTAsSQRTCRMTA) 

GO  TO  5 

4 RMTA«RADTRK<NST) 

5 CONTINUE 
NN«10HOBS  CK 

IF  (OEBUG.GT.O.)  PRINT  12*  NNvAZTAtELTA.RNGTA* ASPTAtFLOAT (K) .FLOAT 
ICKD.RHTA 
L»1 

IF  (AZTA.6T.90.)  L«2 

IF  (RNGTA.GT.RMTA.O.ELTA.LT.ELMINCNSAM.L) .O.ELTA.GT.ELHAXCNST) ) RE 
ITURN 

IF  (IRCNST). GT.O.)  GO  TO  11 

C CHECK  IF  PENETRATOR  IS  IN  RADAR  COVERAGE  OURRING  FLIGHT  TIME 
0ELT».25*(TI-TFIRE) 


51 


NIN>0 

00  10  I>1*3 

TCK»TFIRE*I*OELT 

K»L1 

IF  (NSYS0P.E0*3«A.TCK.LT.TSTART(2).A.K.e0.2)  K=1 
K1»K 

IF  (IR(NST).GT.O)  K1>K1«2 

CALL  COORD  (NST.XSITEINSAM) tOFSETt TCKt AZCtELC«RNGC«ASPC*K) 

IF  (NINTRPRINSTI.GT.O)  GO  TO  6 

APEN>TRP lASPC  •SI6TH(l,Kl)fSI6(lfKl)«NSIG(Kl)) 

RMTC« ( APEW/RTRK ( NST ) ) •• ( .25 ) 

IF  (NINTRPR(NST).LT.O)  GO  TO  7 
GO  TO  8 

6 RMTC*RA0TRK(NST) 

7 IF  (ASPC.GT.ASPHAX(NST) .O.ASPC.LT.ASPMlN(NST>}  60  TO  9 

8 CONTINUE 
L»1 

IF  (AZC.GT.90.)  L=2 

IF  (RN6C.LE.RMTC»A.ELC.GE.ELMIN(NSAM*L).A.ELC.LE.ELNAX(NST))  GO  TO 
1 10 

9 CONTINUE 
NIN«NIN«1 

10  CONTINUE 

IF  (NIN.GE.2)  RETURN 

11  CONTINUE 

C SAM  CAN  GET  OFF  A SHOT 
ISH0T=1 

NN-IOHGOOO  SHOT 
IF  (OEBUO.OT.O.)  PRINT  12*  NN 
RETURN 
C 

12  FORMAT  (1X*A10*12F10.A) 

END 
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SUBROUTINE  COORO  (N«XSlTE«OFSETtTC*AZ*EL«RNGEt ASPf K) 

C 

C 

C THIS  ROUTINE  MAPS  TRAJECTORIES  INTO  3 DIMENSIONAL  SPHERICAL  EARTH 
C COORDINATE  SYSTEM  WITH  SAM  AT  (0t0*8495)KM 

C CALCULATES  AZtEL*RNG*ASP  OF  PENETRATOR  WRT  SITE  FOR  SPHEERICAL  EAR 
C TM 

C 

c 

COMMON  /TRAJ/  NXYZT(2) ,XYZ( 1000 *3*2 ) tT (1000t2> «XSTART*TSTART (2>  vTS 
1T0P(2>  tTINdOfZ)  •NTIN*XECM(20«2)  «NXECM«THIN ( 1 Of  2>  •OXSTART*NXSTARTt 
2XSAMEP*XSAMEW*NSYSOPfTSTOPPfTSTARTPfTARGETYtDXYZOTaOOOf3t2»  f IS(2) 
COMMON  / ASAM/  HRAO ( 1 0 ) f RTRK ( 1 0 ) t ELM AX (l0)*TlNlT(10)f T INTER ( 1 0 ) t NS < 
110)«NSS(10lfTRELOAO(10)  f AVVEL(10)f  ASPMiNdOlfASPMAXdC)  fAZMAX(lO)  • 
2RNG(20fl0)fELR(20fl0)»NRNGd0)  fFUS(20f  10)  fFLF (20f  10)  fNFUSdO)«IRd 
30)tTISHd0>fECMEd0)fALTNINd0)fALTMAXd0)fSIGTH(20f4)fSIG<20f4)fN 
4SlG(4)fRLOCKdO)»XMISL(20«10)*TMISL(20flO)fNXMlSLdO)fRAOTRK(10)fR 
SADLOCKdO) 

COMMON  /FARM/  DFdO)  fOEBUGf FLTWTH(2f 2)  f RLETHdO)  f DEBUGl  t TITLE  (6)  fP 
lKSSdO*2)fCORWOTHfCLUTTER*T£RRANEdOf3f  10)  fNTERAClO)  fNTERfTFRACdO 
2)  tANTHdO)  ,AREAfOEBUG2«OFSETOfRELEASE»VELPENfOLEVdO)  fNOLEVf  AVERAG 
3EfNirLl  (6)  •NirL2(6)  *SYMETRYf TEGRESSt FEGRESSt NITL3 (6)  tNINTRPRdO)  tN 
4ASPfFASTRUN 
COMMON  /INTER/  L 

DIMENSION  US(3)f  UP(3)«  USP(3)*  XYZTC(3)f  UPl(3)f  VEL(3) 

C COORDINATE  SYSTEM  ORIGIN  AT  CENTER  OF  SPHERE 
C X ALONG  CORROOOR 

C Z VERTICAL  THROUGH  SAM  SITE  ^ 

C Y ORTHOGONAL  / 

ASP«0. 

ER>6378.165 
RTD-57, 29578 
ANGal80. 

IF  (AREA.LT.O.)  ANG=OFsETD 
UPHl)aCOS(ANG/RTO) 

UPl (2)«SIN(AN6/RT0) 

USMaUS<3)>ER-»HRA0(N)  / 

US<l)eUS(2)=0. 

CALL  THPl  (TCfXYZTCfK) 

XYZTCM=ER^XYZTC(3) 

SX-XYZTC 1 1 )-XSITE>XSTART 
UP  < 1 )aXYZTCM»SlN(SX/XYZTCM) 

SY»XYZTC(2)*0FSET 
UP (2) =XYZTCM*SIN (SY/XYZTCM) 

UP(3)sSQRT  (XYZTCM*«2-UPd)**2-UP(2)**2) 

RNGE-0. 

DO  1 I«lt3 
USPd)>UPd)-USd) 

1 RN6E«RNGE«USPd)*«2 
RNGEaSQRT(RNGE) 

EL«ASIN(USP(3)/RNGE)*RTD 
AZ«0f 

AZLEN»SQRT<USP(1 )*«2^USP<2)**2) 

IF  fAZLEN.LE.O.)  GO  TO  2 
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00T« (UPl ( 1 > •USP ( 1 ) ♦UPl ( 2 ) *USP ( 2 ) ) / AZLEN 

AZ«ACOS(OOT)*RTO 

AZbAMO0(AZ«720.«360«) 

2 CONTINUE 
ASP»0. 

IF  (NASP.GT.O)  RETURN 

IF  (ASPMAXCN>-ASPPIN(N).GE.180..A.NINTRPR(N),GT.0)  RETURN 

IF  (TC.GT.TCL-l.KH  60  TO  A 

SI6N«1. 

00  3 I«l»3 

3 UPin)>XYZ(LtItK) 

60  TO  6 

A SI6N«-1. 

00  5 Ib1*3 

5 UPHI)»XY2(L-1*I*K) 

6 CONTINUE 
UPIM»ER*UP1 <3) 

SXsUPl(l)-XSITE«XSTART 
UPl  n)«uPiM«siN(Sx/upm) 

SY-UPl I2>^0FSET 
UP1(2»»UP1M*SIN(SY/UP1m> 

UP 1 ( 3 ) *SORT  < UP 1 H**2-UP 1 < 1) **2-UP 1 ( 2 ) •*2 ) 

OOT»VELM»0. 

00  7 1=1.3 

VEL ( I ) *SI6N* (UPl ( I ) -UP ( I ) ) 

VELM«VELM^VEL  < I ) **2 

7 OOTaOOT*VEL(I)*USP(I) 

VELNaSQRTCVELM) 

ASP«1 80 .-ACOS ( OOT/ ( VELM*RNGE ) ) ♦RTO 

RETURN 

ENO 
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FUNCTION  TRP  (XltX*Y»N) 

COMMON  /INTER/  II 
DIMENSION  X<N)«  Y(N> 

IF  (N.LE.2)  60  TO  2 
DO  1 I«2tN 

IF  (Xl.LE.XdM  60  TO  1 

1 CONTINUE 

2 IsN 

3 CONTINUE 

TRP«Y(I-l)^<Xl-X(I-l))*<YII)-Y(I-l)  )/(X(I)-xn-in 

Il*I 

RETURN 

END 
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SUBPOUTINe  TRPl  (Xl«ZtK) 

COMMON  /TPAJ/  NXYZT (2) «XYZ(1000t?*2) t T(1000«2) tXSTART«TSTART (2) ,TS 
1TOP(2>«TIN(10«2) «NTINtXECM(20«2) «NXECM*THIN ( 1 0*2) «OXSTART*NXSTAPT* 
2XSAMEP*XSAMEW*NSYSOP*TSTOPP*TSTARTP*TARgETY*OXYZOT( 1000*3*2) «IS(2) 
COMMON  /INTER/  I 
DIMENSION  Z(3) 

I N«NXYZT<K) 

L»IS<K) 

IF  (X1.LT.T(L-1*K))  60  TO  2 
00  1 IsL*N 

IF  (X1.LE.T(I*K))  60  TO  4 

1 CONTINUE 
I«N 

60  TO  4 

2 DO  3 M»2*L 
I»L-M*2 

IF  (X1.GE.T(I*1*K))  60  TO  4 

3 CONTINUE 
1»2 

4 CONTINUE 
OX«Xl-T(I-l*K) 

DO  S J*l*3 

5 ZIJ)>XYZII-1*J*K)^0X*0XYZDT(I*J*K) 

IS(K)sI 

RETURN 

END 
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FUNCTION  SGN  (X) 
IF  (X)  1«2*3 
SGN«-I. 

RETURN 

SGNaO. 

RETURN 

SGNsl. 

RETURN 

END 
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SUBROUTINE  TIMEIN  <ItOT*N) 


THIS  ROUTINE  DETERMINES  ROUGH  ESTIMATES  OF  THE  INTERVALS  THAT  A 
TRAJECTORY  IS  IN  RANGE  OF  A SAM  SITE 


COMMON  /TRAJ/  NXYZT (2) .XYZ < 1000t3*2) »T ( 1000»2) *XSTART t TSTART (2) .TS 
1T0P(2>  tTINaOfZ)  *NTIN»XECM(20»2)  »NXECM»THIN(  I0«2>  fOXSTARTtNASTARTf 
2XSANEP*XSAMEW*NSVSOP,TSTOPP.TSTARTP*TAR6ETYtDXVZOT(lOOO*3.2).IS(2) 
COMMON  /ISAM/  NTYPE (100) tNSITE ( 100 > «XSITE ( 100 ) vPUP ( 100 1 tELMINi 100* 
12)»NTOTS*SITWDTH(100t3)»AVSHOT(100f3)«TERFR(100)*SITRAO(100)*ARSIT 
2(100) 

COMMON  /ASAM/  HRAO ( 10) »RTRK ( 10) tELMAX ( 10) «TINIT ( 1 0) *TINTER ( 10) •NS( 
110)«NSS(10)  »TRElOAO(10)  »AVVEL(10)  »ASPMIN(10)»ASPMAX(10)  tAZMAXdO)  • 
2RNG(20*10)«ELR(20*10)fNRN6(10)«FUS(20f 10)*ELF(20*10) «NFUS( 10 ) t IR ( 1 
30) «TISH(10)«ECHE(10) • AlTMIN( 10) • ALTMAX ( 10) *SIGTH (20«4) vSIG (20*4) «N 
4S1G(4) tRLOCK ( 10) «XM1SL (20*10) yTMISL (20« \0) tNXMISL ( 10) «RAOTRK ( 10 ) tR 
SAOLOCKdO) 

COMMON  /PARK/  OF (10) »0EB00*FLTMTH(2.2) *RLETM (10) fDEBUGl .TITLE (6) »P 
lKSSd0.2)*CORMOTH.CLUTTER»TERRANEd0.3.10)*NTERA(10)  .NTER.TFRAC (10 
2)*ANTH(10) »AREA«OEBUG2,OFSETD.RELEASE.VELPEN«OLEV(10>«NOLEV.AVERAG 
3E.NITL1(6) •NITL2(6) .SYMETRY.TE6RESS.FEGRESS.NITL3(6) .NINTRPR(IO) .N 
4ASP.FASTRUN 

DIMENSION  VECI(2)f  THMM(2) 

DIMENSION  IN(2).  SVX(2).  SVY(2) 

DIMENSION  XY2T0(3) 

EQUIVALENCE  (XYZTO ( 1 ) .XTO) . (XYZT0(2) .YTO) » (XYZTO(3) .ZTO) 
0PR-57.29S7 

XMlNaXSITE(I)-RLETH(N)-XSTART 
IF  (AREA*LE.O.)  XHINxSlTRAO( I)>RLETH(N) 

XMAX«XMIN*2.*RLETH(N) 

ARSITd)siaO. 

IF  (XMIN.GE.O..A.AREA.LE.O.)  ARSIT d ) =ASlN (RLETH(N) /SITRADd ) ) •DPR 
NTIN*0 

FLTV(TH(1.1>=-1.E*10 

FLTVTHd.2)=-l.E*10 

FLTMTH(2.2)al.E*10 

FLTWTH(2.1)al.E.10 

J»1 

TO*TSTARTP-DT 
KaQ 
INL«0 
1 K»K.l 

IF  (K.GT.2)  K«1 

IF  (K.EO.l)  TO«TO*OT 

IF  (K.EO.l .A.NSYSOP.EQ. 2)  60  TO  S 

IF  (K.EQ.2.A.NSYSOP.EQ.1)  60  TO  S 

IF  (NSYS0P.NE.2.A.K.EQ.2.A.(T0.LT.TSTART(2)*0.T0.GT.TST0P(2) ))  GO 
ITO  5 

IF  (T0.GT.TST0PP..5*DT)  GO  TO  9 
IN(K)aO 

CALL  TRPl  (TO.XYZTO.K) 

IF  (AREA.6T.0.)  GO  TO  3 


1 

j 
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IF  (TO.EQ.TSTART)  80  TO  2 

SVY(KI»YTO 

SVX(K>»XTOT 

2 CONTINUE 
XTOT»XTO*XSTART 

XT0»S0RTI<XT0*XSTART)**2*YT0**2) 

IF  (T0»GT.TSTART>  GO  TO  3 

SVXtK)aXTOT 

SVY(K)«YTO 

3 CONTINUE 

IF  (XT0.«E.XM1N«A.XT0.LE*XMAX.A.ZT0.GE.ALTMIN(N) .A.ZTO.LE.ALTMAX (N 
1)>  IN(K)3l 

IF  (iNfKi.GT.O)  FLTWTH<ltK)*AMAXl(FLTNTH(l*K».YTO) 

IF  tlNIKI.GT.O)  FLTWTM(2tK)»AHINl (FLTWTH(2»K).YT0> 

IF  (AREA.6T.O.O.IN<K)*EO*O.O.INL.6T.O)  GO  TO  A 
VMAG»SQRT (SVX (K ) **2*SVY (K ) **2) 

VECI(l)sSVX(K)/VHA6 

VECII2)=SVV(K)/VMA0 

TMETAI*ATAN2 ( veci <2) ♦ VECI ( 1 1 ) *OPR 

THMM ( 1 > *THMM (2 ) »THETAI 

INL*1 

GO  TO  6 

4 IF  (INL.LE.O.O.AREA.GT.Q.)  GO  TO  6 
VM«XT0T*«2^YT0**2 

IF  (VM.LE.O.I  go  TO  6 

VX»XTOr/VM 

VYsYTO/VM 

00Ta-VECI(2>*VX*VECI(l)*VY 
DANG«90.-AC0S(00T)*0PR 
TMET  A I »TMET  A 1 ♦DANG 
THMNnisAMINI  (THMM 1 1 ) .THETAI ) 

THMMI2)3AMAX1<THMH(2) *THETAI) 

VECItUsVX 
VECIC2l»VY 
GO  TO  6 

5 IN(K)sO 

6 CONTINUE 

IF  (K.E^.l)  GO  TO  1 
GO  TO  (7»8)*  J 

7 IF  (lN(ll.Ea.O.A.IN(2>.EQ.O)  GO  TO  1 
J*2 

NTIN>NTIN^1 

TlN(NTINfl)»AMAXUTO-OT*TSTARTP) 

INL*1 
GO  TO  1 

8 IF  nNn>.GT.O.O.IN(2).GT.O)  GO  TO  1 
J»1 

TIN(NTlNf2)«T0 

INL«0 

THIN (NT IN* 1 1 «THMN ( 1 ) -APSIT ( I ) 

THIN(NTIN«2)aTHMM(2)«APSIT(I) 

GO  TO  1 

9 CONTINUE 

IF  (NTIN.LE.O.O.INL.LE.O)  GO  TO  10 
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TIN<NTIN«2)«TST0PP 

THIN (NTIN* 1 ) sTHHM  < 1 ) -ARSIT ( 1 ) 

THINfNTIN«2)>THMMI2)«ARSIT(I) 

10  CONTINUE 

IF  (NTIN.LE.O.O.AREA.GT.O.)  60  TO  12 
00  11  J>1*NT1N 

IF  (XMIN.OE.O..A.THIN(J«2>-THIN(J«1).LT.360.I  60  TO  11 
TNIN(J«1)30. 

THIN(J*2)>360. 

11  CONTINUE 

IF  (AREA.LE.0..A.NTIN.GT.1)  CALL  COLAPS 

12  CONTINUE 

IF  (OEBU6UGT.O..O.OEBUG2.6T.O.)  PRINT  13«  NTiNt  (TINIK*!)  •TIN(K*2) 
lfTHIN(K«l)«THlN|K«2)tKslfNTIN) 

RETURN 

13  FORMAT  CeH  TlN/OUTt IlOt (4F10.2) ) 

END 
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on  nn  n 


SUBROUTINE  COLAPS 


THIS  ROUTINE  COMBINES  OVERLAPPING  ANGULAR  COVERAGE  INTERVALS 


COMMON  /TRAJ/  NXYZT (2) .XYZ ( 1000«3f 2) tT ( 1 000*2) • XSTARTt TSTART (2) t TS 
1T0P(2) tTINn0*2)«NTIN*XECM(20*2>*NXECH*THINI10*2) tOXSTART*NXSTAPTf 
2XSAMEP*XSAMEW*NSYS0P»TST0PP*TSTARTP«TARGETY*0XYZ0T( 1000*3*2) » IS (2) 
NIbNTIN-1 
N«NTIN 
00  e i3i*Ni 
IP»I^1 
00  7 Js2*N 
K«N*IP-J 

TI1«AMOO(ThIN(I»1)*720.*360«) 

Tl2aAMO0aHlN(I*2)*720.*360.) 

TK1>AMOO(THIN(K*1)*720.*360.) 

TK2aAMOOfTHlN(K*2)*720**360.) 

TIM»<TI2^TI1)*.5 

TI0«(TI2-TI1)*.S 

IF  (TIO.GE.O.)  GO  TO  1 

TIO«180.*TIO 

TIM«AMOO(TIH*900**360.) 

1 TKM»(T»C2*TK1)*.5 
TK0»<TK2-TK1)*.5 

IF  (TKD.OE.O.)  GO  TO  2 

TKO*180.*TK0 

TKMbAMOO ( TKM  *900 • * 360 • ) 

2 CONTINUE 

AMsAMAX1(TKM*TIM)-AMIN1 (TKM*TIM) 

IF  (AM.6T.180*>  AMS360.-AM 
IF  (AN.GT.TK0«TI0)  go  TO  7 
TlN(I*l)«AMINl(TIN(I*l)*TIN(Kf 1) ) 

TlN(I*2)«AMAXi(TIN<I»2)*TIN(K*2)) 

IF  (TKM. 6E. TIN)  GO  TO  3 

S*TKM 

TKMaTIM 

TIH»S 

S-TKD 

TKOsTIO 

TID«S 

3 CONTINUE 

IF  (TKM-TIM.GT.180.)  TIM«TIM*360. 

TMIN(I*1)«AMIN1 (TKM.TKO*TIH-TIO) 

THIN(I*2)>AMAX1 (TKM*TKO*TIM«T10) 

IF  (THIN(I*2)-TH1N(I*1) .LT.360.)  60  TO  A 
THIN(I*l)sO. 

THIN(I*2)>360. 

4 CONTINUE 
NsN-1 
N1>N1-1 

IF  (K.EO.N*!)  GO  TO  6 
L1«K*1 
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L2aN«l 

DO  5 L«L1»L2 

THIN(L-1»1)sTHIn«L.1) 

THlNIL-lt2)sTHlN(Lt2) 

5 CONTINUE 

6 CONTINUE 

7 CONTINUE 
6 CONTINUE 

NTIN«N 

RETURN 

END 
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APPENDIX  B 

SAMPLE  PROBLEM  FOR  SURVIVE 

A.  INPUT  CARD  LISTING 


63 


CARO  IMAGES  FOR  SAMPLE  PROBLEM 
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SO.  100. 

(EOR  CARO) 

PKSS  SINGLE  SHOT  KILL  PROBABILITY  FOR  PENETRATOR 


CARO  IMAGES  FOR  SAMPLE  PROBLEM 
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lEOR  CARO) 

EGRESS  DEFINE  EGRESS  TIME  AND  SHOT  FRACTION 


CARO  IMAGES  FOR  SAMPLE  PROBLEM 
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title 

0*0000 

sample  problem  for  survive  nooel 

OLEV 

0*0000 

defence  level  edits 

.250E*00 

.SOOE'OO 

.TS0E*O0 

.lOOE'Ol 

SAM 

1*0000 

define  parameters  for  type  i sam 

NS 

NSS 

IR 

HRAO  rtrk  hlock  ELMAX 

ALTMIN 

3 

2 

1 

.500E-02  0.  0.  .900E.02  . 

300E-0I 

TINIT 

TISH 

T INTER 

treloao  aspmin  aspmax  azmax 

ECME 

.AOOE*01 

*100E«01 

.100E*02 

.SO0E*O3  .9S0E>02  .I81E«03  .I0IE«O3  . 

lOOE.OO 

RAOTRK 

RAOtOCK 

.300E*02 

*2S0E«02 

MISLXT 

1*0000 

MISSILE  TINE  TO  DISTANCE  PROFILE 

XMISL 

THlSL 

0. 

0* 

.100E«03 

*100E*03 

RNG 

1*0000 

LETHxL  ENVELOPE 

ELR 

RNG 

0. 

•100E*02 

.A50E*02 

*aooF«oi 

.700E«02 

«700E«01 

.900E*02 

*600E«01 

FUSE 

1*0000 

DEAD  ZONE 

ELF 

FUS 

a. 

*300E«00 

.900E«02 

*300E*00 

clutter 

*1000 

ground  clutter  angle 

TERRAIN 

0*0000 

define  TERRAIN  MASK  ANGLE  DISTRIBUTIONS 

TFRAC 

ANTH 

.S00E*02 

*500E-02 

AV  elf 

AV  ELR 

PROS 

.IS0E*Q0 

*2soe«oi 

.SOOE^OO 

.2S0E*a0 

*2t>OE*01 

.soaE«*« 

XECN 

0*0000 

ECK  INTERVALS 

0. 

*200E*02 

.500E*t2 

*100E«03 

PKSS 

0*0000 

single  shot  kill  probability  for  PENETRATOR 

.20 

0*00  0< 

lOO  0*00  0*00  0.00  0*00  0*00  0*00 

0.00 

pxss 

1*0000 

SlNGtC  SHOT  KILL  PR01A61L1TV  FOR  WEAPON 

.10 

0.00 

AREA 

1*0000 

select  rectangular  coordinates 

AVERAGE 

1*0000 

SELECT  EXPECTED  VALUE  METHOD 

CORWDTH 

20*0000 

CORRIDOR  HIOTH 

FASTRUN 

3.0000 

TIME  STEP  FACTOR 

OF 

0*0000 

INDIVIDUAL  DEFENSE  LEVEL 

1.00 

0*00  0. 

.00  0.00  0.00  0.00  0.00  0.00  0.00 

0.00 

stmetry 

1*0000 

select  HOMOGENEOUS  BOUNDARY  CONDITION 

NOSIG 

0*0000 

SELECT  MAX  RANGE  CRITERIA  FOR  SENSOR 

1 

0 

0 

0 0 0 0 0 0 

0 

XYZT 

•*1000 

AIRCRAFT  INGRESS.  IBO  DEGREE  TURN.  EGRESS 

X 

T 

z 

T 
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tiesinTS  AvCHAceo  oven  all  orrsETS  apm  icabll 

(START  > 20.00 

SURVIVAL  PROSABIIIIT  VIIKOUT  ECM  SURVIVAL  PROBABIlITT  WITH  ECM 

DEFENCE  LEVEL  A/C  TO  RELEASE  A/C  COMPLETE  HEAPON/REL  A/C  TO  RELEASE  A/C  COMPLETE  WEAPON/REL  A/C»W  W/0  ECM 

.250  .860  .791  .968  .988  .939  .968  .956 
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